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Abstract. The "external" or "bulk" motion of extended bodies is studied 
in general relativity. Compact material objects of essentially arbitrary shape, 
spin, internal composition, and velocity are allowed as long as there is no 
direct (non-gravitational) contact with other sources of stress-energy. Physically 
reasonable linear and angular momenta are proposed for such bodies and exact 
equations describing their evolution are derived. Changes in the momenta depend 
on a certain "effective metric" that is closely related to a non-perturbative 
generalization of the Detweiler- Whiting R-field originally introduced in the self- 
force literature. If the effective metric inside a self-gravitating body can be 
adequately approximated by an appropriate power series, the instantaneous 
gravitational force and torque exerted on it is shown to be identical to the force and 
torque exerted on an appropriate test body moving in the effective metric. This 
result holds to all multipole orders. The only instantaneous effect of a body's self- 
field is to finitely renormalize the "bare" multipole moments of its stress-energy 
tensor. The MiSaTaQuWa expression for the gravitational self-force is recovered 
as a simple application. A gravitational self-torque is obtained as well. Lastly, it is 
shown that the effective metric in which objects appear to move is approximately a 
solution to the vacuum Einstein equation if the physical metric is an approximate 
solution to Einstein's equation linearized about a vacuum background. 



PACS numbers: 04.20.Cv, 04.25.-g, 04.40-b, 45.20.-d 



1. Introduction 

Newtonian celestial mechanics typically describes the motion of widely-separated 
masses using two types of parameters (see, e.g., [D [I]). These concern either the 
behavior of each body as a whole - the "external" or "bulk" parameters - or the 
details of their internal dynamics. Examples of external parameters are the center 
of mass positions and spin angular momenta of the various masses. The internal 
variables include, e.g., the density and velocity distributions inside each body. In 
typical applications, there is very little coupling between the internal and external 
parameters. As a consequence, one can often compute the center of mass positions 
of each extended body in an iV-body system as though that system were composed 
of point particles described only by their positions and masses. Furthermore, the 
spin angular momentum of each body in such a system can usually be taken to remain 
constant (and does not affect the center of mass motion). Both of these statements are, 
of course, approximate. A more accurate description requires introducing additional 
parameters such as quadrupole moments. These depend on the internal dynamics, 
but in a relatively mild way that often lends itself to simple phenomenological models. 
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The external variables decouple from the internal ones in Newtonian gravity 
largely because the net force and torque exerted by a body's self-field always vanishes. 
There are no self-forces or self-torques in this theorjjj] The instantaneous evolution 
equations for an object's linear momentum, center of mass position, and spin do not 
explicitly involve its self-field. These quantities are affected by the self-field through 
its action on other objects, although this is an effect that takes time to accumulate. 

A priori, it is not clear that similar statements can be made for matter interacting 
with relativistic fields. Such fields carry energy and momentum, so self-forces 
arise generically. This does not, however, preclude an internal-external split of the 
dynamics. The usefulness of such a split does not require that self-forces vanish 
entirely, but only that they do not depend in any essential way on the details of a 
body's internal structure. 

To illustrate this point, consider the motion of a small electric charge in 
approximate internal equilibrium moving non-relativistically in flat spacetime. It has 
long been known that under suitable conditions, the center of mass acceleration a(s) 
of such a charge at time s very nearly satisfief[|] 

ma = Fcxt + 1 'J™^- (1) 

6 as 

Here, Fext is an externally-imposed force, q is the object's total charge and m its 
(bare) mass. The last two terms on the right-hand side of this equation arise from 
interactions with the body's own electromagnetic field. The first of these is "simple" 
in that it depends only on bulk parameters - namely q and a - already required to 
describe the motion of a charged test particle. 

dm, by contrast, has a very different character. In an appropriate approximation, 
it is the self-energy of the charge distribution as it would typically be defined [5]. 
Denoting the electric charge density by pe, 

Smis) . 1 / d^X / d3x' . (2) 

It is clear that Sm depends on the body's internal structure in a nontrivial way. Despite 
this, may be rewritten in the form 

„ 2 oda 

ma = Foxt + o9 3-' 3) 
3 as 

where rh := m + dm is interpreted as a renormalized or effective mass. The same 
assumptions leading to the derivation of ([!]) can also be used to show that dm/ds = 0. 
For a well-behaved extended object, Sm is finite. The mass renormalization effect is 
therefore finite as well. 

Even though the self-force is significant in this example (and depends on nontrivial 
details of the body's internal structure), the final equation of motion involves only 

t The self-force is defined liere as the net force exerted by the self-field. The self-field is, in turn, 
defined in a standard way. See Sect. |2.1| below. Note that this definition of self-force is not the same 
as the perturbative one used in, e.g., [3]! 

§ This equation has been established as a valid approximation only for the acceleration of a physical 
charge (see, e.g., [HIS]). This does not mean that a trajectory with an acceleration satisfying ([TJ for 
all time is guaranteed to stay near the physical trajectory. Many such motions violate the conditions 
under which the equation was derived (even on short timescales), and must therefore be discarded. 
Additionally, there may be neglected terms which lead to qualitatively different behavior over long 
times. Better-behaved equations arise by "reducing order" [J] |6], which changes Q only by an 
amount comparable to the error terms that are already present. This leads to an equation often 
attributed to Landau and Lifshitz in the relativistic case [7]. 
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parameters of the same type as those aheady needed to describe the motion of 
a charged test particle. To the extent that ([ij can be trusted, this means that 
the external variables largely decouple from the internal ones in electromagnetism. 
The center of mass acceleration of an appropriate extended self-interacting charge 
distribution is the same as the acceleration of a monopole test charge moving in 
an effective electric field given by the external one plus |gda/ds (at the particle's 
location). This effective field may be shown to arise naturally as a certain solution to 
the vacuum Maxwell equations [3 |H1 [5] . 

This result can be generalized considerably. Essentially all restrictions regarding 
the charge's size, internal dynamics, and speed may be removed. For almost any 
bounded self-interacting charge-current distribution in flat spacetime, physically 
reasonable linear and angular momenta may be defined that evolve as though they were 
the momenta of an extended test charge (or a pointlike test charge "with structure" ) 
moving in a certain effective electromagnetic field [5]. This effective field satisfies 
the vacuum Maxwell equations near the charge. All effects of the self-force and self- 
torque can be non-perturbatively absorbed into the definitions of the momenta and the 
effective field. Whether or not the internal structure is "effaced" from the external laws 
of motion therefore reduces to a question regarding the nature of the effective field. 
In all but the most extreme systems, the effective field may be shown to depend only 
on bulk parameters like the total charge. Very similar results also hold in generic (but 
fixed) curved spacetimes. The only qualitative change that occurs when introducing 
spacetime curvature is that the quadrupole and higher multipole of a charge's stress- 
energy tensor are renormalized along with its momenta |10j . Analogous statements 
are known for matter interacting with linear scalar fields as well (TUl E] . 

Results of this type greatly expand the scope of - and provide a basis for 
- what has been referred to as the Detweiler- Whiting axiom [HI [H]. It is well- 
known that point particles are incompatible with, e.g., the standard formulation 
of Maxwell electrodynamics (and with general relativity [13]). Despite this, "point 
particle methods" can still be used if additional axioms are introduced into the theory. 
Suppose, for example, that a certain portion of the self-field associated with a pointlike 
electric charge is assumed not to affect its motion. Detweiler and Whiting considered 
this possibility with an ignorable field constructed using a certain symmetric Green 
function [9]. Subtracting this field from the physical one leaves a result which is 
easily calculated and well-behaved. It also satisfies the vacuum Maxwell equations at 
the location of the particle. Substituting this difference field into the Lorentz force 
equation produces the standard Dewitt-Brehme result [T^ [Hj for the motion of a 
self-interacting charged particle in curved spacetime. Similar subtractions were also 
used to efficiently reproduce equations of motion that had previously been derived 
for self-interacting scalar charges as well as uncharged masses in linearized general 
relativity. 

The results of [51 [TT] show that this ability to ignore what is referred to as 
the Detweiler- Whiting S-field is not merely a computational shortcut allowing the 
use of point particle methods in cases where actual point particles cannot exist. A 
very general type of "Detweiler- Whiting axiom" may be rigorously derived from first 
principles for a large class of extended scalar and electromagnetic charge distributions 
moving in fixed spacetimes. This paper uses similar methods to treat the gravitational 
problem. Specifically, it investigates whether the bulk dynamics of an uncharged mass 
in general relativity can be reduced to test body motion in an effective metric (in 
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a nontrivial sensejii]). Related questions have been studied in various contexts using 
the post-Newtonian approximation [I] [121 HEl HZ] ; where they are often referred to as 
"efFacement principles" or demonstrations of the strong equivalence principle. 

The work presented here is motivated more by the types of systems commonly 
encountered in discussions of the gravitational self- force. These discussions typically 
allow the body of interest to move at relativistic speeds in a strongly curved 
background spacetime, but restrict it to be small compared to all scales associated 
with that background. One also assumes that the internal structure of the body does 
not vary too rapidly. Under these conditions - made precise in, e.g., [T51 [12] - an 
equation of motion may be derived that does not depend on any details of the body's 
internal structure. At lowest order, it is just the geodesic equation associated with 
the background spacetime. The next approximation introduces forces due to both 
gravitational self-interaction and spin. The latter effect is the Papapetrou force long 
known to act on spinning test particles [501 [S] ■ The self- force component is typically 
referred to as the MiSaTaQuWa force after the authors who originally obtained it: 
Mino, Sasaki, Tanaka, Quinn, and Wald [331 [3S]. Neglecting the Papapetrou term, 
the motion is most naturally viewed as geodesic with respect to a certain effective 
metric satisfying the linearized vacuum Einstein equation [5]. 

We show that this is a special case of a much more general result. Certain 
definitions of linear and angular momentum are proposed for extended compact matter 
distributions in general relativity. It is assumed that there is no stress-energy near 
the object of interest other than its own (except perhaps dark energy equivalent to 
a cosmological constant). An effective metric gab is then defined based on a non- 
perturbative generalization of the Detweiler- Whiting decomposition of the physical 
metric gab- There is a sense in which the force and torque depend only on gab and 
the details of the body's stress-energy tensor. If gat varies sufficiently slowly that it 
can be expanded in a Taylor series about an appropriate point inside the body (in 
a Riemann normal coordinate system constructed using gab), the instantaneous force 
and torque are shown to be identical to those of an appropriate test body moving in 
the effective metric. A similar result also holds for a certain definition of the center of 
mass. This means that equations known to hold for test bodies (possibly with higher 
multipole moments) also hold for masses with significant self-interaction. 

As a simple application, note that the simplest test bodies move on geodesies. The 
simplest self-interacting bodies therefore move on geodesies of the effective metric. The 
MiSaTaQuWa expression for the gravitational self-force follows easily with some minor 
additional assumptions. Similarly, the simplest equations for a body's spin evolution 
are those of parallel transport. Appropriate self-gravitating masses therefore parallel- 
propagate their spins in the effective metric. Corrections to these statements arising 
from higher order multipole moments are easily added when appropriate. 

The assumptions adopted here are different from those found in other treatments 
of the gravitational self-force. Most importantly, the approaches of, e.g., [THl [TO] are 
intrinsically perturbative. They work in an intermediate "buffer" region outside of 
the body of interest and assume that the metric there is a small perturbation off of 
some vacuum background. Such methods can be applied even to the motion of black 

II Suppose that it is known, for example, that the acceleration of an extended nonrelativistic electric 
charge satisfies ma = qE + f . If g 7^ 0, this is trivially equivalent to the motion of a pointlike test 
body in the field E -|- g~^f. In general relativity, equations with the form mDu'^(s)/ds = can 
always be rewritten as geodesic equations associated with some connection. It is only in special cases, 
however, that such identifications are useful. 
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holes, which Ues beyond the scope of the formahsm developed here. We require that 
a body be described by a well-behaved stress-energy tensor. Despite this restriction, 
there are considerable benefits to our assumptions. They allow the analysis of objects 
that may be highly distorted and dynamical: Dixon's multipole expansions [53] for 
the motion of extended test masses are generalized to all orders. Explicit formulae for 
the momenta are also provided in terms of the body's internal structure. There are 
aesthetic advantages as well. In regimes where they overlap, the method presented here 
requires far less computation than others in the literature. It also provides significantly 
more physical insight. 



Layout of the paper 

The main results established in this work are obtained using only a modicum of 
computation. Despite this, a number of concepts and techniques are employed that 
are not in common use. While not new [51 [TOl HH US] , the relevant ideas are reviewed 
in Sect. [2] by applying them to problems with which the reader might be more 
familiar. Sect. |2.1| starts by discussing in detail the motion of self-interacting masses in 
Newtonian gravity. While the conclusions of this section arc completely standard, the 
formulation used to obtain them has several unusual features. It does not, for example, 
rely on any choice of coordinates. It also treats a body's linear and angular momenta 
as different components of a single scalar functional on the space of Euclidean Killing 
fields. The lack of self-forces and self-torques is shown to follow from the symmetries 
of a particular Green function used to define what is meant by the term "self-field." 

Similar techniques are used in Sect. |2.2[ to discuss the motion of fully relativistic 
extended test masses in curved spacetimes. This is a review of appropriate aspects of 
Dixon's work on the subject [H [231 [H] as reformulated in PF. Linear and angular 
momenta are shown to arise as two components of a scalar functional that now takes 
as input certain "generalized Killing fields" when ordinary Killing fields do not exist. 
These vector fields are defined in detail in the appendix. Using them, multipole 
expansions for the force and torque are established when appropriate. A center of mass 
is also defined, and it is pointed out that the hidden momentum is generically nonzero 
(i.e., the linear momentum is not parallel to the center of mass velocity). Familiarity 
with the discussion of Sect. [2] is essential for understanding the remainder of this 
paper, where similar techniques are used to analyze the motion of self-gravitating 
masses in general relativity. 

The main results of this paper are contained in Sect. [3| After discussing the 
problem of defining an effective metric and a self-field abstractly in Sect. |3.1[ a 



specific definition for the self-field is given in Sect. 3.2 Exact expressions for the force 



and torque arc then derived in Sect. 3.3 Multipole expansions of these equations are 
performed in Sect. |3.4| A brief discussion of the multipole moments appearing in the 
resulting series is given in Sect. [3.5| Finally, a center of mass is defined in Sect. |3.6| 
and related to the linear momentum. 

Sect. |4| applies these results to the motion of a small body. The monopole-dipole 
approximation is discussed in various ways. We then specialize to general relativity 
linearized off of a vacuum background and derive an equation for the center of mass 
position that includes the MiSaTaQuWa "self-force." A similar result is also obtained 
for the spin evolution (including a "self-torque"). 

The sign conventions used here are those of Wald [17]. Metrics therefore have 
signature +2 and the Riemann tensor satisfies 2V[aVb]CJc = Rabc^^d for any 1-form 
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uja- The Ricci tensor is given by Rab — Racb'^- Multiple metrics are discussed in this 
paper, so indices are not raised and lowered unless indicated otherwise. In almost all 
cases, factors of the appropriate metric are displayed explicitly. There are three main 
metrics that appear: gab denotes the full physical metric, gab a certain effective metric, 
and gab a background metric. Derivative operators and curvature tensors associated 
with the latter two geometries are distinguished with a hat or bar as appropriate. 
Non-geometric quantities (like momenta) with hats typically denote renormalized or 
effective versions of their plainer counterparts. Abstract indices are written using 
letters from beginning of the Latin alphabet, while Greek indices represent spacetime 
coordinate components. The letters . . . represent spatial coordinate components. 
Units are used where G = c = 1. 

2. Motion in simple cases 

The main goal of this paper is to describe, in some sense, the large-scale or bulk 
motion of extended masses in general relativity. This is done by analyzing quantities 
that may be interpreted as a body's net linear and angular momenta (as well as the 
closely related notion of its center of mass) . 

The type of momentum considered here is similar to the one developed by 
Dixon [31 [231 [551. Mathematically, Dixon's momenta are tensor fields defined non- 
perturbatively along a preferred worldline in the physical spacetime. They take as 
input this worldline and a timelike vector field prescribed along it. The linear or 
angular momentum of an extended body is then computed by integrating its stress- 
energy tensor over a spacelike hypersurface in a particular way. The evolution of these 
quantities is strongly constrained by stress-energy conservation. 

The only significant restriction to the use of Dixon's momenta is that an object's 
stress-energy tensor be bounded in spatial directions. This bound is also required 
not to be "extremely large" in a particular sense [24l [28] . Limitations on the metric 
are minimal. Despite this, most applications (e.g., [5S1[3D]) have been restricted to 
the test body regime where the body of interest is not allowed to backreact onto the 
geometry. While Dixon's momenta retain a number of interesting properties in a more 
a general context (2 [231 [HI [S2 > other characteristics are less satisfactory. For example, 
it has been shown that even in flat spacetime electromagnetism, the momenta do not 
behave as simply as might have been expected once electromagnetic self-interaction 
is taken into account |33| . This problem can be eliminated with a relatively simple 
modification |S]. 

Similar changes are proposed here in order to obtain physically reasonable 
momenta that obey simple evolution equations in the presence of significant 
gravitational self-interaction (but without electromagnetic or other long-range non- 
gravitational fields). The basic strategy is to first postulate "bare" momenta. These 
agree with Dixon's definitions in the test mass regime, but differ in general. The 
important point is that the evolution equations for the bare momenta include total 
time derivatives of certain terms involving parts of the self-field. These derivatives 
are easily eliminated by redefining the momenta. The resulting variables obey simple 
evolution equations in a wide variety of contexts. 

This section sets the foundation for deriving these results by reviewing the relevant 
techniques in simpler systems. We start by discussing the motion of self-interacting 
extended masses in Newtonian gravity. This is carried out from a somewhat unusual 
point of view introduced in [11] . Similar techniques are then used to analyze relativistic 
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motion in a curved spacetime, but without self-interaction. The resuhing definitions 
and conclusions are equivalent to Dixon's [H |Ml US] . Many aspects of the formalisms 
discussed in this section carry through almost without modification to cases involving 
self-interaction in general relativity. 



2.1. Self-interacting masses in Newtonian gravity 

As a first step to understanding the motion of extended bodies in general relativity, 
consider the motion of a freely-falling extended mass in Newtonian gravity. Such 
a mass may be modelled as having a time-dependent configuration on a three- 
dimensional Euclidean space gab)- We assume that this configuration has nonzero 
volume and may be entirely contained in a compact region C at time s (i.e., 
the body is extended but with a finite size). Also suppose that the body of interest is 
composed of a single material with mass density p(x, s) and 3-velocity v°'{x, s). Given 
any x G M and s g M, p{x, s) > 0. While it is possible to relax this requirement 
considerably, we assume for simplicity that p{x,s) and v°'{x,s) are smooth in both 
of their arguments. Lastly, suppose that there exists an open neighborhood of Eg 
containing no matter other than the body of interest. This ensures that there is no 
direct contact with other objects. 

In general, the density and velocity distributions are constrained by local mass 
and momentum conservation. These laws have the explicit forms 

^+^a{pvn^O, (4) 

OS 

and 

^{pv'^) + Vbipv^'v" + E''^) = -A9"'Vb0, (5) 
OS 

where Vq denotes the Levi-Civita connection associated with the Euclidean 3-metric 
Qab- S"'' = ^^°'''\x,s) represents the body's stress tensor, (j){x,s) the gravitational 
potential, and g°'^{x) the inverse of gab{x). Besides these equations, the body of 
interest is also assumed to be a source for the gravitational field. The potential 
therefore satisfies 

5'^''VaVfc0 - 47rp (6) 

throughout E^. 

The simplest consequence of these equations is that the total mass m cannot 
change. Using dV^ to denote the natural (three-dimensional) volume element 
associated with gab, let 

m:= I p{x,s)dV. (7) 

It immediately follows from Q that m is independent of s. 

Eq. ([5]) constrains the evolution of the body's total linear momentum and 
angular momentum Sa. The linear momentum is typically defined by integrating 
the components pv^ in a Cartesian coordinate system X'^{x). A similar integral also 
exists for the angular momentum. It is important for later generalizations to avoid 
any coordinate choices such as these and instead define the momenta geometrically. 
This can be accomplished by recalling that global linear momentum conservation 
is associated with the translational invariance of Euclidean space. Similarly, global 
angular momentum conservation is related to the rotational invariance of Euclidean 
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space. Translations and rotations together comprise the continuous isometries of 
{■M.,gab)- Generators of these isometries are KiUing vectors. 
Given any Euchdean KiUing field and time s, consider 

P^{s) ^ pix, sy'{x, s)gabix)eix)dV. (8) 

This is a linear functional on the six-dimensional space of Killing fields. It may 
be viewed as returning the component of momentum "conjugate" to If, say, 
C° = d/dX^ is a particular Killing field associated with translations in the X^- 
direction, Pc^{s) is equal to the Euclidean component pi{s) = gabP'^C^ of the body's 
linear momentum as it would ordinarily be defined. Similarly, use of a purely rotational 
Killing field in ([S]) returns a component of the body's angular momentum. 

In general, is equal to a sum of linear and angular momentum components. 
The precise form of this sum may be established by recalling that any Killing field is 
fixed everywhere by specifying it and its first derivative at a single point [27] . Choosing 
a (possibly time-dependent) origin 7s G A^, 

r(x)^{r(%),v6r(7.)}- (9) 

This correspondence may be observed explicitly by noting that in Cartesian 
coordinates X^{x), all Euclidean Killing fields have the form 

Cix) = Cils) + [XHx) ~ x^(7,)]a,e(7s). (10) 

dls) may be chosen arbitrarily in this equation, while 9jC(7s) must be an 
antisymmetric matrix. Note that Ci^) is linear in the "data" {Cils) , djS,'' (js)} ■ This 
is a generic feature of the correspondence ^ , and is unrelated to working in Euclidean 
space. 

Now, it is clear from ([s]) that is linear in ^."'{x). It follows that P^ may always be 
written as a linear combination of ^°(7s) and V6^°(7s). The appropriate coefficients 
are essentially the linear and angular momenta as they would typically be defined. If 
Sa := gab^^ let p° and S"'' = satisfy 

P^is) = p'^(7„ s)S,(7,) + ^S'^\rs,s)Va^b{ls). (11) 

All Killing fields may be generated by varying Sa(7s) amongst all possible 1-forms 
and VaSf, = W[a'^b]{js) amongst all possible 2-forms. Knowledge of P^ for all is 
therefore equivalent to knowledge of and 5'°''. The angular momentum 1-form Sa 
may be extracted from 5''^'' via 

Sa leabcS"^, (12) 

where eabc denotes the natural volume element associated with gab- Both Sa and iS*"^ 
cont ain all angular momentum information in the three dimensions considered here. 



Eq. (12 1 may therefore be inverted, giving S — e S, 



Eqs. ([8]) and (11 1 provide a coordinate-invariant definition for a body's linear and 
angular momenta. Mathematically, p°(7s,s) and 5'°'' (7^,5) are tensors at the point 
7s. In more elementary presentations of Newtonian mechanics, 7^ corresponds to the 
preferred point required to define the angular momentum. It is often taken to coincide 
with the object's center of mass, although this choice is not essential. 

Later sections in this paper makes extensive use of functionals like P^. For a 
relativistic object moving in curved spacetime, the quantity appearing in the 
integrand of ([s]) is translated into an obvious analog involving the body's stress-energy 
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tensor. More difficult to generalize is wliicli must be ciiosen from a suitable space of 
"generalized Killing fields" when ordinary Killing fields do not exist. This is, however, 
a surmountable problem. The relation ( |ll| ) between P^, p", and S'^^ does not change 
at all for relativistic motion in curved spacetimes. 

Even in the Newtonian context, there are advantages to working with rather 
than p° and 5*"^. Most importantly, this functional allows the linear and angular 
momenta to be manipulated simultaneously merely by performing operations on 
scalars. also provides a clear relation between symmetries and conservation laws. 

It is now possible to discuss how a Newtonian body's momenta vary over time. 
These changes can be extracted from changes in P^. Using Q together with ([s]) shows 
that 

^^P^{s) = - £ p{x, s)£^^{x, s)dV, (13) 



where C^(f> denotes the Lie derivative of (f> with respect to S^"-. Differentiating ( 11 ) also 
shows that 

where 7" := d7"/ds and we have used the fact that second derivatives of Killing fields 
vanish in fiat space. Equating the right-hand side of this equation with the right-hand 



side of ( |13[ ) produces evolution equations for both and 5'°''. It is useful to define a 
force F°- and torque iV^ = N^"-''^ such that 

^P^{s)^F^Ea + lN^''V,E,, (15) 
ds 2 

Then 

-^=F'^, (16) 
ds 

= 2p[''7,^l + iV"^ (17) 

ds 

where D/ds denotes the covariant path derivative associat ed with gab- Just as 



knowledge of P^ is equivalent to knowledge of p° and S , (15) provides a one-to- 
one correspondence between dP^/ds and P" and N"'''. The (possibly unfamiliar) term 
involving ^^"7^' in the evolution equation for the angular momentum measures the 
degree to which and 7° fail to be collinear. This term vanishes here if 7^ is chosen 
to coincide with the body's center of mass. In the relativistic context, it rarely vanishes 
exactly. 



Now note that the evolution equation ( 13 ) for P^ is linear in cj). It therefore 
makes sense to discuss the force and torque exerted by particular components of the 
potential. Consider, in particular, the effect of the self-field (j)s- This is defined using 
a symmetric Green function Gs(x,x') — Gs{x',x) that satisfies 

g^'VaybGsix, x') = -Andix, x') (18) 

and vanishes when its arguments are infinitely separated. In Cartesian coordinates 
X^{x), it is explicitly 

Gs(.,.') = |x(,)_'x(.or ^'^^ 
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The self-field (/)s is now defined by 



p{x',s)Gs{x,x')dV'. (20) 



It is clear from ( 18 1 that 

.g^'V^Vfc^s = 47rp (21) 
in Eg. Combining this with ([6| shows that the difference field 

(^:=0-0S (22) 
satisfies the vacuum equation 

g^'Va^b^ = (23) 

in Ss- 

Inserting (20) and (22) into (13) and commuting integrals shows that 

^ = - ^ dVpC^^ dvj^ dV'pp'C^Gs, (24) 

where p' := p{x',s). This equation involves the Lie derivative C^Gs{x, x') of a two- 
point scalar Gs{x,x'). Unless otherwise noted, Lie derivatives of objects depending 
on multiple points are defined in this paper to act on each of those points individually. 
For example, 

C^Gs{x,x') = C{x)^aGs{x,x') + C' {x')Va'Gs{x,x'). (25) 

The integrals involving this term in ( |24[ ) determine the force and torque exerted by 
the self- field (i.e., the self- force and self-torque). These are very simple to evaluate. 
Gg is invariant with respect to all translations and rotations, so 

C^Gs = (26) 

for all Killing fields . All self- forces and self-torques therefore vanish. The momenta 
satisfy 

lp, = -l^dVpC,^. (27) 

Instantaneously, this is the same as the equation satisfied by a test body with density 
p immersed in the (vacuum) potential (j). 

Elementary discussions of Newtonian gravity commonly ascribe vanishing self- 
forces and self-torques to Newton's third law. It is instructive to note that this concept 
is equivalent to the symmetry (26 1 of the Green function used to define the self-field. 
To see this, consider two small volumes dV and dV. In Cartesian coordinates, the 
gravitational force exerted on matter in dV by matter in dV is reasonably interpreted 
to refer to 

p{x, s)p{x', s)d,Gs{x, x')dVdV'. (28) 

Now consider only the first coordinate component of the force exerted by dV on dV' . 
Adding to this the first coordinate component of the force exerted on dV^ by dV' 
results in 

p{x, s)p{x\ s)CqGs (x, x')dVdV' = 0, (29) 

where C° = d/dX^ is a translational Killing vector. This argument may obviously be 
repeated for any translational KiUing field C,°- (i.e., for any satisfying Vf,C° = 0). It 
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follows that the force on dV due to dV is equal and opposite to the force on dV due 



to dV. Considering translational Killing fields in (26) therefore implies the weak form 
of Newton's third law. That Gs is also invariant under rotations implies the strong 
form of Newton's third law: Forces that dV and dV exert on each other point along 
the line connecting them. 



Returning to the main development, (27) provides an exact expression for the 
force and torque exerted on an extended mass in Newtonian gravity. It is not, however, 
particularly useful in this form. It is important to take into account that in many 
practical scenarios, the effective potential </> (or "external field" in this context) varies 
slowly inside Sj. One might therefore expect that (j) could be adequately approximated 
inside the body using only the first few terms in a Taylor expansion. Integrating each 
term of such a ser ies recovers standard multipole expansions for the force and torque. 



Noting from ( 23 1 that (f> is harmonic in , it must also be analytic in this region 
(unlike (j), generically) . The effective field may therefore be expanded in a Taylor series 
about, e.g., 7^ £ Eg. While it is not guaranteed that the resulting series converges 
throughout "Eg, we assume that it does. The Taylor expansion of may be written in a 
coordinate-invariant manner by introducing Synge's function cr{x, x') = a{x' , x). This 
is a two-point scalar equal to one-half of the geodesic distance between its arguments. 
In Cartesian coordinates X*(x), 

a{x,x')^^-\X{x)-X{x')\\ (30) 

Derivatives of a{x,x') may be used as "radial vectors" between x and x'. Holding 
x' fixed, one derivative of a{x,x') at x produces a 1-form at x whose coordinate 
components are the "radial vector" 

a,{x,x')^X,{x)-X,{x'). (31) 

Here, we have used the standard notation ctq := VaC. Using C^a{x,x') = 0, it follows 
that 



°° (—1)" 

n=0 



X g^^'Hls) ■ • •5""'"(7.)>C5V,, • • • Vfe„0(7s) (32) 
for all x' € Es and for all Killing fields 



Inserting ( 32 ) into (|27[) and integrating term by term. 



1 00 1 

^/d^) = - E ;^™'^^-°"(s)^cVa, • • • ^aJhs), (33) 

where m"^ '"" is the "complete" 2"-pole mass moment 

m'^--"(s) := (-W^^n7s)---5'^"''"(7.) 

X / p{x',s)abAjs,x')---<JbAj,,x')dV'. (34) 

It is clear that m'^i """ is symmetric in all of its indices. Many of its components 



do not, however, enter the law of motion (33). To see this, note that for any n > 2, 



multiples of the (inverse) metric symmetrized with any tensor of rank n — 2 may be 
added to m"! """ without affecting dP^/ds. This is a consequence of the fact that 



satisfies the vacuum field equation (23). Using this freedom, each rn"-^'""'"- may be 
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replaced by another tensor that is both symmetric and trace-free. For example, the 
complete quadrupole moment m"'' may be replaced by 



The resulting trace-free moments are the ordinary ones fomid in textbooks. Eq. (33) 
is then equivalent to standard multipole expansions for the force and torque acting on 
an extended mass in Newtonian gravity. 

More explicit equations may be obtained by fixing the point 7^. It is natural to 
do so by choosing this to lie at the body's center of mass. The center of mass is defined 
to be the point about which the first mass moment vanishes: 

m^s) - -g^'ils) I p{x\ s)adls.x')AV' = 0. (36) 



A unique solution to this equation is guaranteed by the assumption that the mass 
density can never be negative and m ^ 0. Choosing 7^ such that = eliminates 
the dipole (n = 1) term in (33). Differentiating (36) with respect to s and using ^ 



also demonstrates that 

p"-™7s- (37) 

This is the ordinary relation between linear momentum and center of mass velocity. 
Note, however, that a similar equation does not remain true in the relativistic case 
(although it is often an excellent approximation). 

Explicit laws of motion may now be written down for the center of mass position 



7s, linear momentum p", and spin Sa- Combining (12|, (14), (33), and (36) shows 
that 



ds 



(7s) (mV;,0(7s) + f]^TO"^-="(s)VfcVc, •••VcJ(7.)) , (38) 



n 

n=2 



1 r. QO 1 



Combining the first of these equations with (37) and dm/ds = immediately provides 
a similar expansion for the acceleration of a body's center of mass in terms of its 
multipole moments and derivatives of 4> evaluated at 7^ . 



In most cases of practical interest, the first few terms in (38) and (|39j) provide 
excellent approximations to the true force and torque. If, e.g., the object of interest 
has a size 0{d) and is separated from other objects by a distance of 0{D), successive 
terms in the multipole expansions tend to differ in size by a factor of at least d/ D. 
A more precise bound may be obtained using standard expressions for the remainder 
term associated with a Taylor series of finite order. Even better estimates can be 
found using Fourier transforms. See, e.g., |34| . 

Note that the gravitational potential (j) entering into the final laws of motion is 
not the one that would be measured using local experiments (which is 0, or really its 
gradient). The body's momenta satisfy evolution equations that are instantaneously 
identical to those of an extended test mass with moments in"'^'"'^" moving in the 
vacuum field </> (^0). 

It is the intent of this paper to demonstrate a similar result for self-gravitating 
masses in general relativity. This is done in two steps. First, Sect. |2.2| considers 
relativistic test masses moving in a prescribed spacetime. The laws of local mass and 
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momentum conservation ^ and ^ are then replaced by conservation of the body's 
stress-energy tensor. No particular relation is, however, assumed to hold between 
the body of interest and the spacetime metric. Once the relevant techniques are 
established, Einstein's equation is "turned on" in the remainder of this paper and 
self-interaction is dealt with directly. 



2.2. Test masses in curved spacetimes 

Consider a relativistic extended body moving in a curved four-dimensional spacetime 
{M,gab)- Associated with this body is a stress-energy tensor — T'"''' (x). Denoting 
its worldtube hy W := supp T"'', spatial slices of W are assumed to be compact and 
to have positive (but finite) 3-volume. As in the Newtonian case discussed above, all 
laws of motion are to be derived from generic local conservation laws. In this context, 
Eqs. Q and ^ are replaced by stress-energy conservation: 

VaT°^ = 0. (40) 

We assume for now that the body of interest is a test mass, meaning that it does 
not affect the spacetime metric gab- There is therefore no replacement for ([6]) in this 
section. 

Now consider the "momentum functional" 

P^(I]) := / g^teT'^dS,, (41) 

where dSc denotes the natural 3-surface element associated with gab- P^C^) takes 
as input a hypersurface S assumed to bisect W and a vector field that is chosen 
later. As with the similar functional Q defined in the Newtonian case, -Pj(S) may be 
viewed as returning the component of momentum conjugate to at a "time" defined 



by E. This interpretation is completely standard if is a Killing vector: Eq. (40) 
then implies that -P^(E) is independent of S (i.e., it is conserved). 

As in the Newtonian case, linear and angular momenta and 5'°'' = S''°''l may 
be defined by demanding that be a linear combination of these two quantities. 



Specifically, it is useful to retain (11). This relation does not, however, make sense 



without being more specific about the types of vector fields that may be used 



in (41). In general, there is no reason to expect that any Killing fields exist (and 
certainly not the 4 -I- 6 = 10 required to define all components of and S""^). Using 



a relation like (11) requires that be chosen from a ten-dimensional vector space 
with the property that each vector is determined throughout a hypersurface given 
knowledge of an arbitrary 1-form 'Rails) and an arbitrary 2-form Va'^b — ^[a^b]ils) 
at one point 7^ € E^. Furthermore, ^°(a;) must be linear in this "initial data." 

There are many spaces of vector fields with these properties. We now specialize 
to specific definitions that recover Dixon's definitions [H [231 for the linear and 
angular momentum of an extended body. Using the terminology of [25' , is assumed 
to be of the fornjj] 

r = (42) 

where 2a is a Killing- type generalized affine coUineation constructed using gab- This 
is defined precisely in the appendix. Following |5l [lOl [11] , we simplify the terminology 



^ The simpler notation = g"'''^b is not used in order to avoid confusion when multiple metrics are 
introduced below. 
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Figure 1. A schematic illustration of the geometry required to define generalized 
Killing fields; The shaded region denotes a portion of the body's worldtube W. 
7s is a point on the timelike worldline F. is a timelike vector at 7s and Ss is 
a spacelike hypersurface formed by the union of all geodesies orthogonal to n" . 



by referring to the (or simply as generalized Killing fields (GKFs) with respect 
to gab- 

Defining GKFs requires fixing not only a metric, but also a timelike worldline 
r = {7s I s G R} and a timelike vector field n° S Tj^J\4 along T. The worldline serves 
as an origin about which to compute multipole moments of T"''. The fix a family 
of spacelike hypersurfaces that provide a time function 9 x i— > s inside the 
body's worldtube W. At fixed s, E^ is defined to be the union of all geodesies that 
pass through 7^ and are orthogonal to n° at that point. These geodesies are not to 
be extended so far that they intersect either with each other (except at 7s) or with 
another hypersurface in the family. It is assumed that the body is sufficiently small 
that such restricted geodesies still form hypersurfaces Eg that foliate W. See Fig. [T] 
for an illustration of the geometry. Under mild assumptions, 7s and can both be 
specified uniquely using center of mass conditions [55] (see also ( 68 ) and ( 69 1 below) . 
For now, however, we continue to describe the general case where they are left free. 

Once gab, T, and n° have been fixed, vector fields = 5"^2b that may be used 
in are to be chosen using the definitions in the appendix. The result is a ten- 
dimensional vector space with a number of characteristics that are very similar to 
those of genuine Killing fields. First among these is the "rigidity property" that has 
already been mentioned: Given any 1-form Sa(7s) and 2-form VaSf, = V[aSf,](7s) at 
a single point 7s G F, a GKF Sa(x) is fixed for all x in the neighborhood W of F 
defined in the appendix. The Sa(a;) are linear in Sa(7s) and VaSb(7s). 

GKFs are also "approximately Killing" near F, meaning that 

V(,Sb)|r = V„V(bS,)|r = 0, (43) 

or equivalent ly, 

C^gablr - Va/:55hc|r = 0. (44) 

The space of generalized Killing fields includes any genuine Killing fields associated 
with gab- If C° satisfies C^gab = everywhere, it is also a generalized Killing field. 
In maximally-symmetric spacetimes, the space of generalized Killing fields coincides 
with the space of genuine Killing fields. The dependence on a preferred worldline and 
foliation disappears in this special case. More generally, the may be interpreted as 
the generators of an "approximate Poincare group" for an observer moving on F (and 
with a preferred time-slicing determined by the Es). 
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Linear and angular momenta p°'{s) and 5"'' — S'[°^l(s) may now be introduced as 
tensor fields along F. Following ( [TT] ), let 

(45) 



Pc(S.) =p"(s)S„(7,) + i5'"'(s)V,Sfc(7,). 



As in the Newtonian case, knowledge of for all possible is equivalent to knowledge 
of p° and S'^^. Here, however, S"'^ is not equivalent to an angular momentum 1- 
form. This second-rank tensor has six independent components in the four spacetime 
dimensions considered here. Roughly speaking, there is a sense in which S"^'' contains 
information associated with both the spin 1-form and the mass dipole moment of 
Newtonian physics. 

It is possible to write Sa(a;) directly in terms of ^aijs) and VaSh(7s) if x G Ss- 
The resulting expressions use Synge's function a{x,x'). Recall that this is a biscalar 
defined to equal one-half of the geodesic distance between its arguments. To borrow 
terminology from optics, it is essentially a characteristic function for the spacetime. 
Many properties of a{x,x') are discussed in [TH [211 ES] ■ As is standard, we denote 
derivatives of a by appending indices: e.g., aaa' ■= VaVa'Cr. Defining °-{x,x') := 
[—craa'{x,x')]~^, and if a{x,x') :— H°- ''{x,x')aab{x,x'), all GKFs may be shown to 
satisfy [25] 

e\x')=g-\j,)[K-\{j,,x')Et{js)-H-''{ls,x')a,{js,x')VtE,{j,)] (46) 
if x' lies within the normal neighborhood of 7^ and x' € Sg- This equation generalizes 



the Euclidean expression (10 1. Inserting it into (41) and (45) provides explicit 



expressions for the momenta as integrals over the body's stress-energy tensor: 



P''is)=g'"'{ls) I ga'b'{x')K'^\{is.x')T''^'{x')dS,,, 



and 



S^\s) = 2 / ga>b'{x')H-'^\-is,x')g''^-{^s)crc{ns,x')T''''\x')AS,,. 



(47) 



(48) 



These momenta coincide with standard textbook definitions in flat spacetime. In 
curved spacetimes, they are the momenta identified by Dixon as being particularly 
useful for the description of objects with conserved stress-energy tensors [5J [211 [2S]- 
Formulae (47) and (48 1 are included here for completeness, but are not needed in any 



arguments below. 

As in Newtonian physics, changes in the momenta may be computed from changes 



in P5. First note that (44) may be used to show that for any GKF and any 7^ G F 



V6VaS,(7,) = Rcab\ls)^d{ls)- (49) 

An equivalent relation holds everywhere for genuine Killing fields [27] . Differentiating 



ab 



( 45 1 while using this identity. 
As is standard, the notation 7° used here denotes the tangent vector to the curve 



F at 7s. Eq. (50 1 provides a recipe for extracting the covariant derivatives Dp°/ds 
and D5'°''/ds from dPc(I]s)/ds. The only difference between this equation and its 
Newtonian equivalent (14) is the presence of the Riemann tensor Rbcd""- This arises 
from ([491 . 
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Everything said thus far has involved only definitions. It is now possible to 
explore the physical properties of the momenta that have just been described. In 



general, P^{Y,s) depends on s. Using (40) and (41 1, the difference in between two 
times s and s' > s is 



T-'C^gabdV, (51) 

/J2(s,s') 

where n{s, s') is defined to be the portion of the body in between the two hypersurfaces 
Eg and E^/ . dV denotes the natural (four-dimensional) volume element associated with 
Qab- Eq. ( |5T| ) may be put into differential form by letting dS :— f^dSa, where f is a 
time evolution vector field for the foliation {Ss}: 

"^-Pd^s) ^ \ f T^'C^gabdS. (52) 



d 



s 



Equating the right-hand side of this equation with the right-hand side of ( 50 ) provides 
evolution equations for p° and S°'^. 



It is clear from (51 1 that is a conserved quantity if is Killing. If dP^/ds ~ 
for all one recovers the Papapetrou equations [201 HI] typically used to model a 
spinning test particle. More generally, changes in P^ measure the deviation from these 
equations. In this formalism, Papapetrou terms in the laws of motion arise purely as 



a kinematic consequence of ( 44 ) and ( 45 ) 



The discussion up to this point has not made any strong assumptions regarding 
the nature of the metric. In particular, self-fields have not been excluded. We now 
assume, however, that in a Riemann normal coordinate system X^(a;) with origin 
7s, the metric components may be accurately expanded throughout Y.^ H W in 
a Taylor series about 7s. In particular, introduce four 1-forms at 7s. These are 
assumed to form an orthonormal tetrad, so 

g''\ls)e^aeb = r)^\ = 3a6(7s), (53) 

where 77^^, = rf"^ — diag(— 1, 1, 1, 1). This tetrad allows the introduction of four 
Riemann normal coordinates X'^(x') associated with the point x': 

X^{x'):^-e^,g'^\jMjs,x'). (54) 



This definition is, in part, motivated by the Euclidean expression (31). Note that 
X^ils) — 0, so 7s is the origin of this coordinate system. 

The metric in Riemann normal coordinates can be viewed as a matrix of scalars 
that depend on the choice of origin 7s and the coordinates X^{x'). This matrix is 
given by [TD] 

gf,dls,X'^{x'j) = v^xV.peXH'''''{ls,x')H'''''{j,,x')ga'b'{x'). (55) 
Taylor expanding these scalar^ in Xf^{x') about X'^ = leads to a general 
expression for the metric that does not make any explicit reference to the e^. Letting 
X°'{'js,x') := — g"^(7s)(T6(7s, a;'), the resulting series is [10] 

ga'b'{x) = Craa'{ls,x)atb'{'-fs,x')g'"'{js)g'"^{'-fs) 
00 ^ 

^ E;^^^'(T-^')---^^"(T-^Wi-/„(7a). (56) 

n=0 

+ In general, g^i, need not be analytic at 7s. We assume that a finite powe r se ries nevertheless 
provides an adequate approximation throughout Eg n W. The Taylor expansion ( |56| should therefore 
be cut off at finite n. We write an infinite upper limit and an exact equality sign here (and in similar 
equations below) for simplicity. 
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The gab,ci---c„ils) appearing in this equation are referred to either as tensor extensions 
of gab or as metric normal tensors. They are derived from the coefficients appearing 
in the Taylor series for g^^^: 

Despite appearances, the gab,ci - c„ do not depend on the choice of e^. The zeroth 
extension is the metric itself and the first extension vanishes. In general, it is clear 
that the n"^ metric normal tensor is symmetric in both its first two and its last n 
indices. It may also be shown that [TU] 

5a(&,ci---c„) = 5(a6,ci--c„_i)c„ = (58) 

for all n > 2. Keeping this restriction on n, all metric normal tensors can be written 
as polynomials in the Riemann tensor. To linear order [24|, 



9ab,ci---Cn 



^^l ^ '^{c3-c„{R\a\ciC2)''9bd) + 0{R^). 



(59) 



This equation is exact for n 

^abc 



2, 3. For higher n, there are additional terms nonlinear 
Using certain details of the GKFs together with (56), one may derive a power 



in Rnh/^ or its derivatives 



series expansion for C^gab [10] 
^iga>b'{x') = 



oo 



-t(- 



n=2 



} a' b' d-i---dn 



abci 



X X-^i (7„ x') • • • X"^" {-fs,x')C^gab,c^-c^ hs). (60) 

The omitted coefficients in this series are known explicitly in terms of tr if x' € Ss, 
which is the only case relevant in this section. More generally, such a series still exists. 



although the coefficients are no longer known exactly. Substituting ( 60 1 into ( 52 ) now 
shows that 



oo 

- y ^i^'- 

n=2 



^''\s)C^gab,c,-cSls). (61) 

■^""^(s) appearing here are interpreted as 



By analogy with (33 1, the coefficients I"^^' 
the 2"-pole moments of at the time s. 

Without loss of generality, the symmetry properties of the metric normal tensors 
allow the to be chosen such that they are separately symmetric in their first 

n and last two indices. They may also be taken to satisfy 

j(ci---c„a)b _ jci (c2---c„ab) _ g {^i) 

A unique formula linking moments with these properties to T"^ may be derived using 
( |60| and ^ TI) (see also gl]). Like ^ and (|48]), the resuh has the form of an 
integral over involving the stress-energy tensor and various bitensors constructed 
from a. It is significantly more complicated than the Newtonian formula (34) for 
m""^'""". This is partially because T"^ has two more indices than p. Much less obvious 
is that the relativistic moments are "reduced" with respect to ( 40 1 . They are adapted 



to describing conserved second-rank symmetric tensors. Knowing all of the ' 
together with p° and S'"'' is equivalent to knowledge of T"'' [Mj . The same statement 
does not remain true if is replaced in all integrals by a second-rank symmetric 
tensor that is not divergence-free. 



Mechanics of extended masses in general relativity 



18 



The given index symmetries imply that I'^^' 
(?i + 3)(n + 2)(n- 1) 



has a total of 



(63) 



algebraically independent components. This far exceeds the number typically ascribed 
to the 2"-pole moment in other formalisms [37]. The reason for this is essentially 
that the I " are "complete" in the sense described in the previous paragraph. If no 
restrictions are placed on gat, nothing further can be said. Recall, however, that traces 
of TTi'^i---"" decouple from the Newtonian equation (33) because <?!> is a vacuum field. 
Similarly, certain components of I " decouple from (61) if gab satisfies the vacuum 
Einstein equation Rat — 0. This may be seen by noting that certain traces of (59) 
vanish in this case. Use of (44) shows that these same traces still vanish if Rat = ^gab 
for any constant A. In most cases where a test body description is appropriate, the 
/■■■ may therefore be replaced in (61) by moments with many fewer components. 
Additional discussion of these points may be found in [TU] , although precise details of 
the reduction process are not known. 

Another important point to note is that the sum in (61) starts at n 
corresponds to quadrupole order. It is a consequence of (44) and (52) that the 
monopole and dipole moments of T"'' 
contribute to dP^/ds. 



2. This 

ffld ([52) 

essentially p'^ and S"''' - do not directly 
These moments do, however, affect Dp°/ds and 'DS°'^/ds via 



the Papapetrou-like terms appearing in ( 50 ) . 
a net torque N"'' = N^''''\s) such that 

Dp" 

^ " 



iRbc. 



a Qhc',d 



pa 



Explicitly, define a net force F°(s) and 

(64) 
(65) 



Comparison with (45) and (61) shows that 



E 



lf^-f-^''{s)Wbgcdj,-fAls), 



(66) 



and 



^"''(^) = 2^ ^^^'"■""'^(^)[5//.c....c„(7s)4" + ^5d/Ac....c„_,(7s)4:]5'"'(7.). (67) 

n=2 

The hope in writing these series is, of course, that adequate approximations may 
be obtained by truncating them at some small maximum n. This can only happen 
if r and {Ss} are chosen appropriately (if it is possible at all for a given system). 
We now fix a particular worldline and foliation that is hopefully "appropriate" in 
this sense. This is done by imposing center of mass conditions as described in, e.g., 
[51[2HI3H]- First recall that is constructed using geodesies that pass through 7^ 
and are orthogonal to at that point. Suppose that F and n° are chosen such that 



5a6(7s)p"(s)5^^(s) = 0. 



(68) 
(69) 



Under mild assumptions, the resulting F and exist, are unique, and are timelike 
pS] , The first of these equations essentially states that the {S^} foliation is the one 
preferred by "zero-momentum observers." Eq. (69) encapsulates the notion that 
a body's center of mass position is the point about which its mass dipole moment 
vanishes in the zero-momentum frame. 
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Unlike in Newtonian gravity, the center of mass velocity 7" of a relativistic mass is 
not necessarily proportional to p° . Relating these two quantities is simpler if the time 
parameter s is chosen such that 5afc(7s)p"('S)7s = where the mass is defined 

by 

Ms) [-5a6(7sK(s)/(s)]i/2. 

This means that in general, 7° does not have unit norm, 
of generality in assuming that gabnin'l = — 1. Hence, 



(70) 

There is, however, no loss 



(71) 



With these conventions, differentiation of (69) may be used to show that the linear 
momentum and center of mass velocity are related via |38j 



mis 



P 



ab 



Nf'ghruDRcdb'gfi] 



m2 + ISb-SdfR.Jgfi 



(72) 



This equation breaks down if + \S^'^S'^^ Rhcd Qji — 0, which may be interpreted 
as a constraint on S°'^/m. Such a restriction is implied by the conditions required for 
the center of mass to exist as a unique timelike worldline. 

Note that the center of mass velocity does not appear on the right-hand side 
of (72). Indeed, the complexity of this equation arises mainly from the nontrivial 



operations required to solve explicitly for 7°. 

evident that Eqs. (64)-(67) and (72) form ordinary differential equations (ODEs 
and 5*°^ with the form 



Displaying ( 72 ) in this way makes it 

for 



7s, P 



Is 



Bp'' 
~d7 



BS 



ab 



ds 



(73) 



The right-hand sides of these equations involve only geometric quantities, 7^, S"''', 
and the higher moments /■". It follows that if the metric is known and the quadrupole 
and higher moments are prescribed functions of s, the motion is uniquely determined 
by specifying initial values for 7^, and S""^. There is, however, no physical reason 
that the higher moments should be treated as given functions of s. Generically, their 
evolution depends on the details of the specific system under consideration. 

To summarize, Dixon's momenta have now been described as they apply to a 
compact body with a conserved stress-energy tensor. Eq. (41) constructs from 
and (45) Unks p"^ and S'"'' to Pj, and the are built from the metric in 

the appendix. Taken together, these relations define Dixon's momenta in a very 
general context. Their changes are described exactly by the integral ( [52| (even if the 
body's presence strongly influences the geometry). As in Newtonian gravity, integral 
expressions for the force and torque are not particularly useful by themselves. Ideally, 
one would like to be able to constrain the motion without detailed knowledge of 
the body's stress-energy tensor. This is easily accomplished in the formalism just 
described if the metric does not vary rapidly inside the body (specifically, on a given 
EsHW^). dP^/ds can then be expanded in the multipole series (61 ). This is, however, a 
reasonable procedure mainly for test masses. In cases involving significant self-gravity, 
gab generically varies rapidly inside the body. In such cases, it cannot be accurately 
approximated using a low-order Taylor series. 

Nevertheless, the remainder of this paper establishes that there is a sense in which 
useful multipole expansions can be performed even in the presence of a significant 
self-field. This is because forces and torques exerted by the (potentially complicated) 
self-field are shown only to make an object's momenta and higher moments appear 



Mechanics of extended masses in general relativity 



20 



slightly shifted from what might have initially been expected. In this sense, the self- 
field may be eliminated from the instantaneous laws of motion by appropriate (and 
physically reasonable) redefinitions. 

Obtaining this result requires defining precisely what "should" be meant by 
"linear and angular momentum" as well as "self-field." The momenta used below 
reduce to the definitions provided in this section in a test mass limit. More than this, 



the form (41) for -Pj(S) is retained exactly as written. The space of vector fields from 
which the are to be drawn does change, however. 



3. Mechanics of self-gravitating extended bodies in curved spacetimes 

The remainder of this paper considers an extended mass with stress-energy tensor 
moving in a spacetime {A4,gab)- Its worldtube is denoted by W. A worldline F is 
taken as given together with a collection of hypersurfaces {Ss|s € K} that foliate W. 
Each slice "Eg OW is assumed to be compact and 7^ := F n (see Fig. [T]). Both F 
and {Ss} may be fixed using center of mass conditions (see Sect. |3.6| below), although 
we do not require this. Unlike in Sect. |2.2[ Einstein's equation 
1 



Rab - ^gabR + ^gab = ^ngacgbdT" 



(74) 



is now required to hold with stress-energy tensor T"'^ at least in a neighborhood of W 
(i.e., the body is assumed not to be in direct contact with any other source of stress- 
energy). The presence of a cosmological constant A in Einstein's equation does not 
significantly change any arguments below, so we allow it to be nonzero. Additionally, 
note that T"-^ is automatically conserved as an integrability condition for ( 74 1 . 



The "bare" momentum associated with a body of the type just described is defined 
by introducing a P^{J^s) computed via (41) as a certain integral of T°-^ over S^. All 
metrics and volume elements in that formula are those associated with the physical 
metric gab- This is to be distinguished from a certain "effective metric" gab (7^ gab) 
defined below, gab plays an analogous role in the laws of motion to the Newtonian 
effective potential (j> introduced in Sect. 2.1 

It will be seen to be most natural to choose the appearing in (41) to be 
generalized Killing fields associated with gab- In particular, let 

r = r'Sfc, (75) 



where g"^^ is the inverse of gab and is a generalized Killing field constructed as 
described in the appendix, but with the gab used there replaced by (jab- The worldline 
F is chosen as a set of origins for the GKFs. This is assumed to be a timelike curve 
with respect to gab- Similarly, the are assumed to be hypersurfaces formed from 
the set of all geodesies passing through 7s € F and orthogonal to a timelike vector 
field n° at that point (notions of geodesic, orthogonality, and timelike all being with 
respect to gab)- 



3.1. Initial considerations 



It follows from (40| and (41) that dP^(I]s)/ds may be computed using (51). In most 



cases of interest, the metric varies rapidly inside the body. It is therefore not useful to 
expand the right-hand side of (51) in a multipole series like (61 1. We instead proceed 



as in the Newtonian analysis of Sect. 2.1 Define a self-field using an appropriate 



Green function and compute the effect of this self-field on dPj(Ss)/ds. The self-field 
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should be chosen such that the force and torque exerted by it are "ignorable" in 
an appropriate sense. Subtracting the self-field away from gab should also leave an 
"effective metric" gab that can be well-approximated by a Taylor series in many cases 
of interest. 

Quite generally, one might suppose that gab can be reconstructed from an effective 
metric gab and a "self-field" H"-'' using local algebraic operations: 

gab = Gabi9,H). (76) 

Assuming a relation of this sort, it follows from ( 51 ) that the force and torque depend 
on an integral involving 

T'^'C^gab - [Aab'^ig, H)Cfgcd + Babcd{g, H)C^H"'] , (77) 
where Aab'"^{g,H) and Babcd{g,H) are determined by Gab{g,H). 

The first term on the right-hand side of (77) has a simple effect on dP^/ds. 
Suppose that the are to be derived from generalized Killing fields associated with 
gab as de scri bed above. Comparison with (51) then shows that the term involving 
Aab'^'^ in ( 77 ) contributes forces and torques identical to those exerted on a test mass 
with a stress-energy tensor 

f''':=Vg/^T^''A,d''\g,H) (78) 

moving in a metric gab- Here, g/g denotes the proportionality factor relating the 
volume elements associated with gab and gab- 

^abcd = y/gjl^abcd- (79) 
In any coordinate system, this is equal to the square root of the ratio of the 
determinants deig^u and detg^i^. The result is, however, independent of any 
coordinate choice. 

If C^cjab can be accurately expanded in a Taylor series on a cross-section n W , 
the first term on the right-hand side of (77) contributes forces and torques that may 
be expanded in a multipole series like (61 ): 

^P^{'^,)^\Y.^{ir-----^^-\s)Cfgab,c,...cM + \ I T'^'BabcdC^H^'dV. (80) 
The 

gab.cf -Cn are tensor extensions of gab as explained in Sect. |2.2[ The multipole 
moments (/') " appearing here are the same as the I used in ( |61[ ) with the 
replacements gab — > gab and T"'' T°'^ . These shifts may be interpreted as (finitely) 
"renormalizing" the body's quadrupole and higher multipole moments. 

It is the goal of this paper to show that there is a sense in which all of dP^/ds 
(or really dP^/ds for an appropriate P^) can be expanded in a multipole series. 
Demonstrating this requires specializing further so that the second term in ( 80 ) may 
be simplified. At this point, gab, H'^^, and Gab{g,H) have not yet been specified. It 
would be ideal if precise definitions could be provided such that 

(i) The self-field H"^'' does not affect the force and torque in any "essential" way. It 
may shift the momenta (as is well-known to occur in electromagnetism) and the 
higher moments, but do nothing else. There should therefore exist some £^{s) 
and some (/")ci - c„a6^^-) ^^^^i that 

I I T'^'BabcdC^H^'dV = I £ ^{ir'---^-'^'{s)C^gab.c.....M - (81) 

•^^'> n=2 ' 

The £({s) appearing here is to be interpreted as a "self-momentum," and must 
depend on properties of the system only in a finite neighborhood of E^ . 



Mechanics of extended masses in general relativity 



22 



(ii) gab satisfies the vacuum Einstein equation Rab — ^gab as closely as possible in 
a neighborhood of W. Following the Newtonian analysis of Sect. |2.1[ one might 
expect that this implies "slow variation" of the effective metric in typical systems. 

Unfortunately, it is not clear how to realize both of these requirements exactly. 
We choose to implement the first precisely as stated. The second is weakened to 
demanding that gab approximately satisfy the vacuum Einstein equation linearized 
about an exact vacuum background gab if dab is sufficiently close (jab (although the 
definition of gab is non-perturbative and does not require a choice of background or 
that gab be "close" to a vacuum solution in any sense). 



3.2. Defining a self- field 

Taking cues from perturbation theory in Lorenz gauge (see, e.g. 



]), choose the 



Gabig, H) appearing in (76 1 by supposing that gab is equal to "a background" cjab phis 
a trace-reversed perturbation H'^^: 



gab — gab + [gacgbd ~ ^gabgcd)H' 



cd 



(82) 



Further suppose that the self-field H is related to the body's stress-energy tensor 
via 



ab 



where Gf" {x,x') 



w 

^(ab)a'b' 



G" aba b' 
S 



ga'c'gb'd'T^''''dV', 



(83) 



G. 



G. 



ab(a'b') . 



is an appropriately-chosen Green function 



associated with metric perturbations about gab- 

Again using Lorenz-gauge perturbation theory as a guide, let Gf satisfy 

□Gf + 2{5^("i?,,/) - J [g'^\R,f - Ag,f) + [rg'^'Rch - A.r'').9d/] }Gf 

= -A,,g^'^'f^'g^\,g'^^a'5{x.x'). (84) 

Here, □ := g'^'^^ d and g°'a'{x,x') may be any two-point tensor that reduces to 
a Kronecker-(5 in the limit that its arguments coincide. For definiteness, we take 
g'^a' to be the parallel-propagator associated with gab. Given any vector at x' , 
g°'a'ix,x')v°' returns a vector at x that is equal to v"" parallel-transported to x using 
gab (along the g-geodesic connecting x and x'). If dV is used to denote the volume 
element associated with gab, the Dirac distribution S{x,x') satisfies 



S{x,x')f{x')dV' = fix) 



(85) 



M 



for any test function f{x) and any x G M. It follows from (83) and (84) that H"^'' 
satisfies the differential equation 



OH^b + 2{g^^'^Rd,f'^ - \ [r\Rdf - Agdf) + {rg'^'Rch - ^'''Ygdf] ]H^f 



ab 



(86) 



If desired, (82) may be used obtain a similar equation for gab as well (written in terms 
of geometric object s as sociated with gab)- 

Note that Eq. ( 84 ) does not define Gg''° completely. There are many solutions 
to that equation (for a fixed gab)- A particular solution may be fixed by making the 
following choices [T?| : 
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(i) Like the Gs used to define the Newtonian self-field in Sect. 2.1 demand that the 
gravitational Green function be symmetric in its arguments: 

Gf''''''ix,x') = G'^'''''\x\x). (87) 

This is essential for most of the self-force to vanish due to cancellations that can 
be interpreted as a manifestation of Newton's third law (as far as it can be said 
to hold in this context). 

(ii) Also demand that Gg'"' {x,x') vanish if its arguments are timelike-separated 
with respect to gat- This ensures that the self- field H"''' defined by (83) does not 
depend on knowledge of the body in the distant past or future. 

Enforcing both of these requirements fixes the Green function completely for a given 
gab- The result is (at least if Rat = A = 0) commonly referred to as the S-type 
Detweiler- Whiting gravitational Green function associated with the metric gab 13 [E] . 

If X and x' are sufficiently close, it is known that Gg^° ^ (x, x') has the Hadamard 
fornix [H] 

Gf = ^[g'^=5^'^5(°'c/)<jAi/2j(CT) - f'^'"^''''e(<7)]. (88) 

Here, S and 8 are the Dirac and Heaviside distributions, respectively. (t{x, x') = 
a{x',x) is Synge's function computed using gab, A(x,x') = A{x',x) is the van Vleck 
determinant, and g°' a the parallel propagator. A is defined to be the unique biscalar 
satisfying 

g'^'aaWblnA^A-g-'aab (89) 

and A{x,x) = 1. This equation can be viewed as an ODE along the g-geodesic 
connecting x and x'. In coordinates, the solution is 

Aixy) = - ''''!^-t^\ (90) 
V-9V-9' 

The tail V"'''"''' (x,x') = ^ {x' , x) appearing in (88 1 is a certain homogeneous 
solution of ( pi] ): 

- \ irf'Rch - A5°').9<i/] V^^'"'"' = 0. (91) 
The appropriate boundary conditions are found by integrating the ODE 

r^Oc^dV'"'"''''' + - 2)l>'^'"''''' = ]^r''^,SJ^(g-Sg^'-g(-' ff)t,A^/^) 

+ [g'^^'Rdcf"^ - \{rg'^Rcp - ^'''Ygdf 

- \g''\Rdf - A5,/)]5'V5''^'h/^AV2. (92) 

along all g-null geodesies emanating from x' . V"-^'^ ^ may therefore be found by solving 
a characteristic initial value problem. 

A Green function satisfying all of the given equations is uniquely determined if 
its arguments are not too widely separated (there should, in particular, exist exactly 

* Ref. [12) only derives this result in cases where Rab = A = 0. The derivation is easily extended to 
the case considered here with no change in the conclusion. 
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one geodesic connecting those arguments). It is important to note that the trace of 
(x, a;') is always proportional to the metric. Contracting (84) with gab, 

1 



aba'b' 

s 



DigabGf^ " ) + (4A - -R) {gabCf^ " ) = -4^^ " 5{x, x') 



(93) 



It follows from this together with the imposition of appropriate boundary conditions 
that 

^a'b' / 



- /^aba'b' 
gab^S 



(94) 



where Gs{x,x') is the S-type Detweiler- Whiting Green function associated with a 
massive, nonminimally-coupled linear scalar field. Gs satisfies 



(4A--i?)Gs 



-At:6{x, x') 



(95) 



and is symmetric: Gs{x,x') = Gsix',x). This Green function also vanishes when its 
arguments are timelike-separated with respect to gab- Explicitly, 



where 



Gs = i[Ai/2<5(a)-fe(a)], 



t7 ^ {'/-aba'b' 

V - -^gabga'b'V 



(96) 



(97) 



Given a physical metric g^b and a body with stress-energy tensor T"*", the effective 
metric gab is taken to be a solution of the simultaneous equations (82), (83), (88 1, (91), 
and (92). This definition is highly implicit, and it is not clear that any solution exists 
(or that it is unique). Despite this, we assume that a unique solution of the given 
system does exist and that it has the same signature as gab- 

One practical method for computing the effective metric is via iteration. As a 
first guess, suppose that gab ~ g'^ab ^'^'^ some metric gab ■ Substituting this zeroth-order 
solution for gab into ( 82 1 and ( 83 ) produces the first-order approximation 



where 



,(1) 



^(i) 

ilab 



jj(Q)ab 



gab Wac ilbd 2^°''' ^'^d. )^ ' 



G. 



{Q)aba'b' .(0) ^(0) rpc' d' 



W 



9a'c'9b'. 



(98) 



(99) 



gal,' may, in turn, be substituted back into (82 1 and (83 1 to find a second-order 
approximation for gab- This process can be repeated as often as required. It is not 
clear that such a procedure converges uniquely for reasonable choices of g^^ij , although 
we assume that it does. 

The definition of the effective metric provided here was chosen essentially as 
the simplest non-perturbative generalization of the Detweiler- Whiting decomposition 
[SI IH] typically presented for solutions of the linearized Einstein equation. It follows 
that gab has particularly nice properties in the linearized regime. Temporarily 
suppose that gab is an approximate solution to Einstein's equation linearized about a 
background metric gab satisfying the vacuum equation Rab = ^(jab- Define a trace- 
reversed metric perturbation = H^"'') such that 

gab = gab + [gacgbd - \gabgcd)H'"^- (lOO) 
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If T-L"^^ is sufficiently smalj^] and gab is an exact solution to Einstein's equation with 
stress-energy tensor T°^, 

an"'' + 2g''^''Rdcf''^n'^f' + {g^^Va ~ 2S'f g''^ f V f)V cV.'"^ = -mnT"'' + 0{H^). (101) 

One may now solve for gah using the iterative method described above. It 
is natural in this case to use g gp as a zeroth-order guess for the effective metric: 
9^ab — 9ab- Combining (98) and (100) then shows that the first-order approximation 
for (jab is given by 

1 

T 



ilab 



9ab + {gac9bd — ^9ab9cd)'HR , 



(102) 



where n"^ n"'' - Combining the first-order analog of (86) with (101), 



(5°'Vd - 2s^;g''^fVf)v,n^^' = oin' 



(103) 



Comparing ( 102 ) and ( 103 ) with ( 100 ) and ( 101 ) shows that g'^^ is approximately a 
solution to the vacuum Einstein equation linearized about gab- It is also expected in 
this linearized regime that 5^^^ is an excellent approximation to gab- It follows that the 
effective metric itself is very nearly a vacuum solution if gab is approximately a solution 
to Einstein's equation linearized about a vacuum background. Note that although 
the definition of the self-field was inspired by expressions for metric perturbations 



in Lorenz gauge, no such gauge choice is required for this conclusion. Eq. ( 103 ) is 



the linearized Einstein equation for trace-reversed metric perturbations in any gauge. 
Furthermore, the equations defining gab are completely non-perturbative. They do 
not depend on any choice of background. 



3.3. Force and torque 

Consider, once again, the momentum functional Pj(Es) associated with an extended 
mass. The (exact) difference between a component of momentum at time s versus the 
same component at time s' > s is given by (51). Using this together with (82) and 



( 83 ) shows that 



P5(E.0-^C(S.) 



T"^''{x)C^gabix)dV 



[ dV [ dV'J'^ix^x'), 



where 



J^^{x,x') := 2r°^r"''''4"^ (gachd - ^9ab9cd)ga'c'gb'd'G 



cdc d 
S 



(104) 



(105) 



r2(s, s') C W represents the portion of the body lying between the hypersurfaces 
and (see Fig. jij), while denotes a "partial Lie derivative" that varies x but 
not x' . As discussed in Sect. 3.1 the first term on the right-hand side 



of (104) 



contributes forces and torques that are easily understood if gab varies slowly inside 

tt Technically, one should consider a smooth one-parameter family of metrics gab{^\^) a-^d the 
associated perturbations 'H"'^{x; X). For fixed A, gab{x;X) is a solution to Einstein's equation with 
stress-energy tens or T °''(x: \). If 3^6(2:; 0) = Qabix), it is clear that ■H°-^ix;0) = T''''{x;0) = 0. An 
exact solution of llOli (omitting the 0{'H?) term) is really a solution for d'H.°'^{x;f)) / dX. Similarly, 
the T"** appearing in that equation should be dT"^ {x; 0) / dX. 
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Figure 2. An illustration of the spacetime regions used in |104[ l and | |106| . The 
shaded volume denotes the portion of the worldtube where information is required 
in order to compute £^{T,s) and f^(Sg/). 



the body. The second term in this equation is more comphcated to understand, and 
may be interpreted as the effect of self-interaction. 

essentiaUy represents the force exerted on matter at x by matter at x' . 
For any such biscalar (even if (105) docs not hold), note that 

[ dV [ dVJ'^{x,x') f '^^ f dV'[J'^{x,x')+J^^{x',x)] 

+ \ I dV ( dV'[F^{x,x')-F^(x\x)]. (106) 
^ Jn Jw\n 

if all integrals can be commuted (see Fig. [2|. The first line of this equation may be 
physically interpreted as averaging the force on matter at x due to matter at x' and 
vice-versa. In this sense, it measures the failure of Newton's 3rd law as discussed in 
Sect. O Using (ISTI), (|9l, and (pSl), 



-\F^{x,x')+F^{x',x)] = T-'T'^ " C^{gac9M9a'c'gb'd'Gs 



cdc! d' 



^gatga'b'Gs). (107) 



The Lie derivative appearing here is the ordinary one acting on both x and x' . If 
is an exact Killing vector associated with gab, Eq. (1071 vanishes exactly. In many 
other cases of interest, it is very small. 

The second line of (106) effectively renormalizes P^. To see this, first note that 
1 



dV 



n{s,s') 



w\n{s,s') 



dV'[J^^{x,x')-T^{x',x)]=£^{j:s)-£d^s'), (108) 



where 



f^(S,) ■.= '^ j^^dV j^dV [F^{x,x') - F^{x' 



(109) 



The (four-dimensional) portion of W in the future of Ss is denoted by E+ in this 
equation, while the portion in its past is denoted by Ej. An explicit formula for £^ is 
easily obtained by combining (105) and (109 1. Using the notation 

1 



hab[R] := ^gacgbd 



\gabgcd) I ga'c'gb'd'T-''''Gt''''''dV' 



(110) 



for any spacetime volume R. 



T'^'C^hablKW^ T'^'C^Kbi^tW). (Ill) 
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Despite appearances, this depends on the behavior of the body only in a finite four- 
dimensional region around n W. Recalling that G^"^ ^ {x, x') = when x and x' 
are timelike-separated with respect to gab, portions of W that are timelike-separated 
from all of Eg n do not contribute to (111 I. In simple cases, £{(Es) depends on 
the body's state for times of order its diameter into the past and future. This is the 
shaded region in Fig. [2j 5^ is interpreted as a "self- momentum" conjugate to 
Further discussion may be found in Sect. |3.5| below. 



An exact equation may now be written down for the force and torque acting on 



an extended body with a conserved stress-energy tensor. Substituting ( 106 )-( 108 1 into 



(104) and taking the limit s' — >■ s shows that 
d 



ds ^ 2 



dST 



ab 



^i9ab 



dV' T 



a'b' 



gac9bdga'c'9b'd 



where 



w 

'cdc' d' 



'abf/a'b'Gs 



(112) 



(113) 



Eq. (112) is analogous to the Newtonian expression (24) with two differences. First, 
the left-hand side involves a time derivative of + £^ rather than of by itself. The 
presence of an additional term here is interpreted as being due to the inertia of the 
body's self-field. That a massive object must carry a field with it as it moves affects 
the momentum that naturally appears in its laws of motion. The other difference 
between the present case and the Newtonian one is that the Lie derivative of the 



Green function appearing in (112) does not necessarily vanish. As explained in Sect. 



|2.1[ that term may be viewed as measuring the violation of Newton's third law. 

As a whole, the component of "effective" force or torque dP{(Es)/ds associated 
with a particular vector field measures the degree to which gab is preserved by 
transformations in the direction To see this 
purely geometric manner from gat 
functional in C^gab- It therefore follows from (112) that dPj(Es)/ds depends on 



note that Gg^° ^ is constructed in a 

\aba'b' 



This means that C^C^ is a (nonlocal) linear 



only via terms linear in C^cjab- This dependence is confined to a finite neighborhood 
of Es n sufficiently large to allow the computation of G'^"' ^ {x,x') for all null or 
spacelike-separat ed p airs {x,x') with x G 'Eg D W and x' g W. 

In general, (112) may be viewed as providing a recipe for computing dP^/ds = 
d{P^ +£^)/ds in terms of C^gab- The earlier Eq. (51) expresses dP^/ds in terms of 
^^Qab- These relations are mathematically equivalent. No aspect of Einstein's equation 
is required to transform one into the other. No properties of the are required eit her 
All that is needed is the decomposition of gab into cjab and H""^ described in Sect. 3.2 



The value of ( 112 ) over (51 ) lies in the physical expectation that C^gab is "better- 
behaved" than C^gab- This expectation derives from the detailed definition of gab, 
and does depend on Einstein's equation. As shown at the end of Sect. |3.2[ gab is 
an approximate solution to the vacuum Einstein equation at least if gab is sufficiently 
close to a vacuum metric. One might therefore expect that in this regime, most of 
the details of a body's internal structure necessarily present in gab do not appear in 
gab- A finite power series expansion for C^jjab is therefore likely to be more accurate 
than one for C^gab (given reasonable choices of ^''). A finite multipole expansion of 



the integral in (112) is then more useful than a similar expansion of (51 ) 
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The construction presented here is by no means claimed to be the only reasonable 
one. There is considerable freedom to change the definition of gab without spoiling 
any of the properties just described. One can, for example, produce definitions 



where the equivalent of (112) includes a three-point interaction term involving 
£^[Gf"^ {x,x')Gg'^ '' {x' ,x")]. An appropriate modification of this sort could be 
extremely valuable. It might, for example, be possible to define an effective field gab 
that is an exact solution to the vacuum Einstein equation (unlike gab) and such that 
an appropriate dP^/ds depends only on L^gab- Whether or not this is possible to 
accomplish in a physically reasonably way is a question that must be left for later 
work. 

The definitions presented in this paper are "minimal" in that they are essentially 
the simplest non-perturbative generalization of the decomposition used in the 
Detweiler- Whiting axiom for the motion of point masses [HI HI] . They also produce 
equations that closely match those that are already known for non-gravitating objects 
with electric or scalar charge [Sl|TTj. It is, unfortunately, untrue that gab is an exact 
solution to the vacuum Einstein equation. This does not mean that the definitions 
used here necessarily fail to be useful for systems where the nonlinearity of Einstein's 
equation cannot be ignored. The interior of, e.g., an isolated star with constant density 
has a metric that is (as expected) quite simple in Fermi normal coordinates based at 
the center of the star Even for a very compact star of this type, the metric 

components are accurately approximated by low-order power series. It appears very 
likely that a similar result holds inside similarly-simple stars in, e.g., binary systems 
(before tidal disruption). These comments apply to the full physical metric gab- The 
effective metric gab is likely to have an even simpler behavior than gat- Precise details 
regarding the properties of the definitions presented here are left for future work. 

3.4- Multipole expansions 



Given T and gab, (112) provides a generic (and exact) prescription for the force 
and torque acting on an extended body. If there exists an exact Killing vector (^^ 
associated with the effective metric gab, 

C^gab = C^Gt"'"' = C^Gs = 0. (114) 

The effective momentum P; associated with is therefore conserved. This is 
analogous to the result that the bare momentum Pq conjugate to is conserved 
if Ci^gab = 0. 

Generically, there are no Killing fields associated with gab or gab- The integral 



appearing in (112) is then rather unwieldy to evaluate directly. It would be 
significantly simpler if C^gab{x) could be expanded in a Taylor series about 7s that 
provides a good approximation for all a; in a sufficiently large (but finite) ball enclosing 
Tis CiW. That such an expansion simplifies the first term on the right-hand side of 



(112) is clear from the discussion in Sect. 2.2 That it also simplifies the second term 
follows from the aforementioned fact that C^G^"^ ^ is a linear functional in C^cjab- 

We now assume that the effective metric can be accurately expanded about 7^ 
in a Taylor series in a Riemann normal coordinate system with origin 7s (where the 



coordinate system is constructed using (jab)- This means that, by analogy with (56) 



ga'b' — CTaa'CTbb'g 9 } , — 
n=0 



ac.bi^ lxf^---Xf-g,dj,...u, (115) 
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denote tensor extensions of gab (see Sect. 2.2) 



where X° := —g°-^(j}j and the gab,c_ 
Although the sum here is written as an infinite one, there is no guarantee that the 
series will converge. We instead assume that the right-hand side is an asymptotic 
series for the left-hand side. There then exists a certain n beyond which the series 
should be truncated in order to obtain an optimal approximation. Despite this, we 
continue to display infinite upper limits and exact equality signs for simplicity. 

Recall that the are taken to be GKFs with respect to gab- The origins used 
to construct these fields lie on the worldlinc F. Using this and (115), it is possible to 



establish the "hatted analog" of ( 60 ) 

oo 

C^ga'b'ix') = ^( 



h' d\---dn 



abci ■■■Cn 



X"' ■■■X''-C^gab,c,-cAls)- 



(116) 



n=2 



The omitted coefficients in this equation are complicated but calculable. Details of 
their properties are discussed in |10j . That the sum cannot include contributions from 
n = 0, 1 follows from gab,c — and 



(117) 



Substituting (116|) into \112\ and integrating term by term shows that there exist 
tensors l^^---'^-n^i 



\s) such that 

-Pe(s.) - - ^ -r--^-'^\s)c,gab,c....,M. 

71=2 



(118) 



Note that this is identical to the test mass expression (61) for dP{/ds with the 
replacements 



9ab 9ab, 



Pi ^ Pi. 



jci---Cnab y jci - ■ -Cnah 



(119) 



This means that if the series expansion ( 115 ) for ijab is valid, the momenta associated 
with a self- gravitating compact body behave as though they were the momenta of a test 
mass with multipole moments J'^^ ^ '^riab j^ioving in a metric gab- AH direct effects of 
the self-field H'^^ have been absorbed into the definitions of and /'^i - '^"'^*'. This is 
made more explicit in Sect. |3.6| below, where laws of motion are obtained for effective 
linear and angular momenta p° and 3°"^- The arguments are essentially identical to 
the ones given in Sect. |2.2|for extended test masses. 



Note that the multipole series (118) is useful only if an adequate approximation 
is obtained by truncating it at some small finite n. Roughly speaking, this occurs 
if all significant length scales associated with gab are much larger than the body's 
diameter. This is likely to be the case if gab is very nearly a solution to Einstein's 
equation linearized about a vacuum background, gab is then very nearly a solution 
to the linearized vacuum Einstein equation. More generally, it is not clear how useful 
the multipole expansion may be. 



3.5. Effective multipole moments 

Before continuing to describe the motion of an extended body, we first note some 



properties of the renormalized multipole moments appearing (118 1. The effective 
monopole and dipole moments are contained in P^ . They differ from the bare moments 



contained in by the £^ defined in ( 111 ). Explicit formulae for the quadrupole and 
higher effective moments I " appearing in (1181 are not derived here. They may, 



however, be determined (at least approximately) using the techniques developed in 
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To] , This requires comparing (112 1 and (118), and is not difficult conceptually. It 



does, however, require a great deal of tedious calculation. 

Regardless of any detailed relations between the moments and T"'', it is clear 



from (1181 and the discussion of metric normal tensors in Sect. 2.2 that the effective 



moments have the same index symmetries as their bare counterparts This means 
that /ci - c„ab separately symmetric in its first n and last two indices. It also satisfies 
(62 1 with the substitution I " — > /■ ■. Note as well that in the regime where gab is 



approximately a vacuum solution to Einstein's equation, certain traces of the multipole 
moments decouple from the laws of motion as discussed in Sect. 



2.2 and in TIP. 



The physical meaning of the self-momentum £^ can be made considerably more 
transparent by specializing to the case of a stationary system that is a solution to 
Einstein's equation linearized off of a Minkowski metric (jab- Assume that T"''' is 
invariant under the action of a particular time-translation vector field r" :— d/dX^ , 
where X° is a globally-inertial time coordinate for gab- The 5° (constructed using ijab) 
are then approximately equal to Minkowski Killing fields Consider, in particular, 
the Killing fields independent of (i.e., exclude boosts). Then, 

1 



T'''hab[W]eAS, 



(120) 



to lowest nontrivial order in the metric perturbation. 

Now specialize further so that the only significant component of T"^ is 
proportional to {d / dX^) ® [d / dX^) . Also take to be a hypersurface of constant X'^ 
that is unbounded in every direction. Suppose as well that the only significant metric 
perturbation is determi ned b y H'^^ (so there is no "external field"). This means that 
^ab _ ^ab _ ^{Q)ab r^^^ effective "energy" is therefore given by 



Pr 



r 



1 



00 



d-^X 



to second order in the metric perturbation. This is equivalent to 



-p, = y[(-g)r"" + <°"]d-^x 



(121) 



(122) 



in the given approximation, where i"'' is the Landau-Lifshitz tensor in the background 

9ab'- 

16^r^ := g.dg^^'VfH-'VhH''' + ^gcd^'V fH^'^^ ^.H^f 



(123) 



Eq. (122) is the usual expression for the energy used in post-Newtonian theory. 



3.6. Center of mass motion 

At this point, effective linear and angular momenta p°'{s) and S'"'' = S'[°''l(s) may be 
introduced as tensor fields on the worldline T used to define the GKFs. By analogy 
with the test mass equation (45), suppose that 

P5(E,) = r mails) + ls-''{s)^a^bhs) (124) 
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for all GKFs Sa(x) defined using the metric gab, the worldline F, and the vector field 
n°. The ap pea ring on the left-hand side of this equation is related to via (75). 
Furthermore, (49 1 holds with the replacements Va Vq and Rabc^ Rabc^- Eq- 



(124) therefore implies that 



ds 



ab 



ds 



2p[''7,^l V.Sfc. (125) 



This is directly analogous to the test mass relation (50 1. Combi ning it with (1181 
yields niultipole expansions for Dp°/ds and DS'^'^/ds identical to (64)-(67) after the 
replacements 



P 

^ab 



p 

■;:,ab 



^ jc\---Cnab y jci - --Cnab 



g- g- 



TD d , D a 

^abc ^ ^abc 



abc 

D/ds 



D/ds, 



9ab,ci---Cn ^ 9ab,ci- 
pa ^ pa ]\jah _ 



(126) 



ab 



We refer to the resulting equations as the "hatted forms" of their counterparts in the 
theory of extended test bodies. 

The worldline and foliation used to construct the GKFs may now be fixed by 
choosing them such that 



p%s) cx <, 

gab{ls)r{s)S'"'{s) = 0. 



(127) 
(128) 



These are obvious generalizations of the center of mass conditions ( 68 1 and ( 69 ) . Unlike 



in that case, however, there exists no proof that (127) and (1281 have well-behaved 
solutions. We assume, however, that they do. 
As in Sect. 



2.2 



it is useful to choose the parameter s such that gabP"'Ts — ~™i 



where the effective mass is defined by 



We also set gabn'!,n\ 



-gabP^'p']"'. 

— 1, so p°- — 



(129) 



m^. The center of mass velocity is then given 
by (72 1 with the replacements (1261 and m rh. Together, the hatted versions of 



(64|-(67) and (72 1 strongly constrain the evolution of the body's linear and angular 



momenta as well as its center of mass. They do not determine it completely. As in the 
test body case, the evolution of the quadrupole and higher moments must be specified 
using other methods. Additionally, the effective metric gab couples to the motion in a 
nontrivial way. This is the main complication in practical computations involving the 
gravitational self-force. 

It can be useful to define a spin 1-form Sa via 



Sa 



leabcanlS^". 



(130) 



The center of mass condition ( 128 1 guara ntee s that all information contained in S""^ 

(131) 



is also contained in Sa- This means that (1301 is invertible 



ibcd ^ Q 



gab ^ ^, 

Note that p'^Sa = 0. The hatted form of (|65|) implies that 



ds 



QabSc 



Bp" 



1 



:^abcdN 



cd 



(132) 
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The Dp'^/ ds appearing on the right-hand side of this equation may be ehminated using 
the hatted form of (64 1. By not doing so, one may interpret the first term in (132 1 



as being responsible for a kind of Thomas precession. It arises from the requirement 
that and Sa remain orthogonal. 

There is nothing that prevents the effective mass m from varying. It immediately 



follows from the definition (1291 that 

dm ^ „ Dp'' 
-T- = -9abn, — . 



(133) 



Substituting the hatted form of ( 64 ) into this equation and simplifying with the hatted 
form of (72) leads to a (large) equation that does not explicitly involve 7". Another 



useful form is 



ds 



= gab -i^F" 



9cd- 



Dp** 
d7 



(134) 



which follows from (128) as well as the hatted versions of (64) and (65). Additional 



manipulations to the right-hand side of this equation may be used to bring it into a 
form involving total s-derivatives and "induction terms" that depend on derivatives of 
the moments in a certain non-rotating reference frame [2l I26j. Regardless, multipole 
expansions for F° and iV"'' start at quadrupole order. Multipole expansions for dm/ds 
therefore start at quadrupole order as well. 



4. Motion of a small mass: monopole and dipole approximations 



As a simple application of the laws of motion just derived, consider the motion of a 
small body around a much larger one. It is possible to adopt precise approximation 
schemes for such problems and proceed rigorously (see, e.g., [TSl [Hj). We instead 
present what is essentially a plausibility argument. This is straightforward to improve, 
although the details are not particularly interesting. 

First consider truncating all expressions for the motion at dipole order. This 
corresponds to ignoring the quadrupole and higher moments Roughly speaking, 
this is a good approximation if all significant length scales associated with gab are 
sufficiently large compared with the body's own size. If gab is nearly a vacuum metric 
(as occurs in the regime of linearized gravity described in Sect. 3.2) there is also a 



sense in which ignoring the higher moments corresponds to assuming that the body 
of interest is nearly spher ical. Regardless of the precise reason for ignoring the higher 
moments, it follows from ( 118 1 that all of the remain constant in this case. Applying 



(125) and the hatted forms of (64) and (65) implies that the force and torque vanish: 
pa _ -^ab _ Combining this result with (134) immediately shows that the effective 
mass rh remains fixed in this approximation. Furthermore, the linear and angular 
momenta evolve via the Papapctrou equations in the effective metric: 
Dp° _ 1 
~d7 ~ 2' 



2p'°7; 



(135) 
(136) 



Using the hatted form of ( 72 1 , the center of mass velocity 7° is seen to be related to 
p"" via 

2 \m^ + \Sb-SdfR,Jgji)' 



P 



(137) 
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Eqs. (135)-(137l form a coupled set of ODEs for p", S"-'', and 7^ . Alternatively, one 
may replace the evolutfon equation for S'^'' with one for Sa using (1301 and (132). 

The laws of motion sir aplif y con sidera bly if the spin can be neglected. Note that 
S"-^ = is one solution to (1361 and (137). It is therefore consistent to consider non- 
spinning bodies in the monopole-dipole approximation adopted in this section. Using 
( 135 ) and ( 137 ) , the center of mass of such a body moves on a geodesic of the effective 

D7s" 



metric: 



ds 



0. 



(138) 



The accuracy of this equation may be estimated by evaluating, e.g., the contribution 
of quadrupole terms to the force and torque using the equations described above. As 
explained in Sect 



3.2 



c/ab is very nearly a vacuum solution to Einstein's equation if 
Qab is itself well-approximated by a solution to Einsteins equation (with stress-energy 
tensor T"'') linearized about an exact vacuum solution. In this context, it appears 
reasonable that all scales associated with gat are sufficiently large compared with those 
of T"'' in, e.g., fairly generic binary systems. Neglect of the higher moments is then 
relatively straightforward to justify. It is less clear what occurs in the highly nonlinear 
regime. It is likely that ( 138 ) remains a good approximation in many cases, although 
it may br eak d own for strongly self-gravitating objects with large inhomogeneities. 

Eq. (138) assumes that the body's spin vanishes. If S"'^/rh is nonzero but still 
small compared to length scales associated with gab, Eqs ( 135 )-( 137 1 may be simplified 
by linearizing in the spin. It then remains true that = rhj'^. Using ( 135 1 and ( 132 ), 
Sa is found to be parallel-propagated along F with respect to gab- 

t>Sa 



0. 

ds 

Substituting ( 131 ) into ( 135[ ) shows that the center of mass velocity satisfies 

67- 1 
ds 2 



-[Kbci e ■> gdhlsls. 



S 



f 



(139) 

(140) 
They 



in this case. It is unclear how consistent (1391 and (140) are over long times, 
do not preserve the constraint p'^Sa — 0, which means that the spin 1-form Sa fails to 
remain equivalent to the spin tensor S"^^. Al tern atively. 5*°^ may be evolved instead of 
Sa using dS'''''/ds = 0. Coupling this with (140) eventually leads to violations in the 
center of mass condition (128 1. Such complications do not arise if Eqs. (135l-(137) 
are retained in full. 

The laws of motion just derived do not explicitly make any assumptions regarding 
the strength of a body's self-gravity (although they are more likely to be accurate when 
the nonlinearity of Einstein's equation can be neglected). They do not require any 
choice of background or gauge. More explicit expressions can be obtained by choosing 
a vacuum background gab and a gauge as well as assuming that gab is "close" to (jab- 

In this context, it is common to compute a body's acceleration with respect t o Qr,h 
rather than gab- This requires writing the derivative operator D/ds appearing in ( 139 ) 
and (1401 in terms of D/ds. For any vector v'^{s) or covector ojai-s), the appropriate 
transformations are 

as as as as 

where 



c:bi>c, 



(141) 



be 



1 



g'"'{2^(^hgc)d-^d9bc)- 



(142) 
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To linear order in the metric perturbation gab 



9ab, 



D7s D7s , 1 



(143) 



ds ds ' 2' 

Note that in general, the background acceleration need not be orthogonal to the 4- 
velocity with respect to gat- This is because the parameter s has been chosen such that 
7" has unit norm with respect to gat (when p"" ~ is negligible). It is more typical 
to normalize the 4- velocity with respect to gat- 
such that u° := (ds/ds)7^ satisfies (jabu'^u^ = - 



Introduce a new time parameter s(s) 
-1. Then, 




TsTs9bc 




(144) 



For sufficiently small spins and weak metric perturbations, the background 



acceleration is obtained by substituting ( 140 1 into this equation 



Tlu" 1 



1 



... ^_ 2-m(r'' + u''u'^)(Vd56c - 2yb9cd)u'u^ + -Rtci'^r'^''^ 9dhU^u' S f . (145) 

The first group of terms on the right-hand side of this expression is what is typically 
referred to as the gravitational self- force in the literature [TH [TSl (THl HH [13]. Note, 
however, that this terminology differs from the non-perturbative notion of self-force 



adopted elsewhere in this paper. The second term on the right-hand side of ( 145 1 is 



an approximation for the ordinary Papapetrou force in the background spacetime. 

There are more "bulk parameters" describing an extended body than only its 
center of mass position. Continuing to define quantities in terms of gab, one might 
be interested in a "background mass" ifi := [—gabP°'p'']^^^- While it has already been 



noted that the m defined by (129) remains constant in the current approximation, the 



same is not true of m. This evolves via 



d_ 
ds 



1 



[gab - gab)u°'u'' 



(146) 



to first order in the metric perturbation. 



The spin evolution equation ( 139 ) may be rewritten in terms of the background 



transformed into ( 145 1 



metric in the same way that Eq. ( 140 ) for the center of mass acceleration is 



ds 



1 



r''{^dgab'2V^agb)d)u''Sc. 



(147) 



Generalizing the terminology typically applied to (145), the right-hand side of (147) 
might be said to represent a "self-torque." This usage is, however, intrinsically 
perturbative. It differs from the notion of self-torque used in other parts of this 
paper. 



Eqs. ( 145 ) and ( 147 ) describe the center of mass and spin evolution of a body with 



small spin and negligible higher moments if gab hes sufficiently close to a background 
metric gab- More can be said by assuming that gab also lies near gab and by adopting 
Lorenz gauge. This gauge choice corresponds to ensuring that the appearing in 



(100) satisfies 



ab 



0. 



Using this in ( 101 ) shows that to linear order in the metric perturbation, 
on"'' + 2g<''Rdcf''^H'^f' = -167^T"^ 



(148) 



(149) 
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The physical metric may be found by specifying data on an initial Cauchy surface 



Esq and using the retarded Green function 0"^^^ {x,x') to evolve this data into the 
future (where Gf^!^ ^ satisfies ( 84 ) with all hats changed to bars) . Then 



n 



ah 



ab 



(150) 



where Hq'' is some homogeneous solution of ( 149 ) and is the portion of W in the 
future of Esq . 

Far outside of W, a body with slow internal dynamics produces a metric 
perturbation nearly indistinguishable from a retarded solution to the linearized 
Einstein equation with a point particle source (see, e.g., [HI UHl IH])- To a first 
approximation, this "effective particle" can be taken to have the worldline F, mass to, 
and no higher moments. This means that 



H°''(a;) ^ 4to 



^rct 



ha'h' 



(2:,7.')ffa'c'(7.')ff6'd'(7.')"''(s')"'''(s')ds' +^^"(0;) (151) 



at large distances. Now, the motion is determined by the difference field 'H'^ appearing 
in ( 102 ). This is a homogeneous solution of the linearized Einstein equation in Lorenz 
gauge, so it may be written in Kirchhoff form: 



'f'9 {G. 



d'f 



Vc'Gf°'^'H'*'^')d'Sh' 



(152) 



Here T is any closed hypersurface enclosing the point x at which the left-hand side is 
evaluated. 

Allowing r to be a very large tube with timelike sides surrounding W, use of 
( 151 ) shows that 



njab 



4to 



r" 



aba b 



ga'c'9b'd'U' u" ds'+Hf, 



(153) 



where C^""'^' 



r^aba b 
^rot 



Q^a b jg typically referred to as the R-type Detweiler- 



Whiting Green function (even though this satisfies a homogeneous wave equation, 
and is therefore not strictly a Green function). The conclusion of this argument is 
that for a sufficiently small particle with slow internal dynamics, the effective metric 
inside W is essentially that of a point particle. A similar statement cannot be made 
for the retarded field. 



The Hadamard form ([88|) for Gg* 
used to compute exp 



aba'b' 



(and the equiv alent for G°g° ^ ) may now be 



icitly. As can be seen from ( 145 ) , we need only the first 
derivative of this field on F. This has already been computed in, e.g., [T^] for the case 
where Rab = and Du°/ds = 0. We shall continue to assume that the acceleration 
is zero, as any terms involving it will be negligibly small. We do, however, generalize 
the Ricci tensor to be Agab- Then 



Vcn^\-/s) = 4to 



H, 



ah 



ab 
5 



(154) 



where 



H, 



ab 



4to lim 



Vcatr''''^^a'Wh'ds'. 



(155) 
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For simplicity, indices in these equations have been raised and lowered with the 
background metric. Also note that the limiting process used to define H°'^c avoids the 
singularity in the retarded Green function. 

Now suppose for simplicity that Wc'Hq'^ is negligible, as can be arranged if linear 



perturbation theory may be trusted sufficiently far in the past. Substituting ( 154 1 



into ( 145 ) then yields 



^ = lir"" + u-u'')u'u'{Hdbc - 2Hbcd) + \Rbci''e'>^''^-gdhu''u'Sf. (156) 

Excluding the spin term, this is the MiSaTaQuWa equation as it is usually written 
(at least if sq -oo) [El [HI UHl ES] . Ref. [ig (which includes the spin term) 



refers to (1561 as the "self-consistent" equation of motion. This is to distinguish 
it from the acceleration of a deviation vector for the center of mass worldline away 
from a background worldline that is geodesic with respect to the background. Such 
constructions have not been used here. 

Lastly, note that an analog of the MiSaTaQuWa equation for the spin evolution 



is easily obtained by substituting (154 1 into (147) 



= -2mu\'Rabc''Sd ~ -u'S^H.ab - 2H(a6)c). (157) 



5. Discussion 



This paper non-perturbatively defines linear and angular momenta and S""^ adapted 
to the study of extended objects in general relativity. An effective metric gab is 
also constructed non-perturbatively from the physical metric gab by generalizing the 
Detweiler- Whiting decomposition originally introduced in [5]. Using stress-energy 
conservation, p° and 5*"^ are shown to evolve via the Papapetrou equations (written 
using gab) plus corrections depending on the degree to which the effective metric fails 
to be symmetric with respect to a certain set of generalized Killing fields. 

Considerable simplifications arise in cases where gab can be accurately expanded 
in a power series throughout the body of interest (in the sense described in Sects. 



2.2 and 3.4). When this occurs, corrections to the Papapetrou equations may be 
expanded in series depending on a body's quadrupole and higher multipole moments. 
These series are formally identical to the multipole expansions provided by Dixon 
O [24l [38] for extended test masses: There is a sense in which self-gravitating masses 
move like test bodies. The metric in which such a fictitious test body appears to fall 
is gab rather than gab- 

The conclusions of this paper may be viewed as providing a justification for the 
gravitational Detweiler- Whiting axiom [3] that a "point mass" moves on a geodesic in 
an effective metric produced by subtracting a certain "S-field" from the physical metric 
(if "point mass" is replaced by "mass with small but finite size" ) . The validity of this 
type of statement has also been extended considerably. It applies to all multipole 
orders and also to the evolution of a body's angular momentum. This joins similar 
results that have recently been established for objects with scalar or electromagnetic 
charge moving in fixed background spacetimes O [H] . 

Future work 



The momenta proposed in this paper are not intended to be the final word on the 
subject. The multipole expansions for the force and torque derived in Sect. 3.4 
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require as their main assumption that gab vary slowly inside objects of interest. That 
this is likely to occur in "typical" systems is motivated in Sect. |3.2| by noting that in 
linearized gravity, the effective metric is very nearly a solution to the vacuum Einstein 
equation. This suggests (but does not strictly imply) that details of a body's internal 
structure necessarily present in gab disappear in (jab- The vacuum condition is also 
useful in that it decouples many components of the "complete" multipole moments 
from the laws of motion. 

It is less clear what occurs in the fully nonlinear regime. In general, gab is not an 
exact solution to the vacuum Einstein equation. Despite this, the multipole expansion 
(118) remains valid in this case as long as gab does not vary too rapidly. Indeed, 
a "second-order self-force" could be derived from this formalism for systems where 
the effective metric varies slowly. Slow variation appears likely inside strongly self- 
gravitating objects with a nearly uniform internal structure, although its precise range 
of validity is not clear. Additionally, many more components of an object's multipole 
moments are required to describe the motion if Rab 7^ ^gab- 
It would be ideal if the definitions provided in this paper were modified so that 
extended objects could be shown to fall like test bodies moving in an effective metric 
that satisfies the vacuum Einstein equation exactly. As has been emphasized at various 
points, there is considerable freedom in the techniques developed here. Different 
choices could have been made in this work without spoiling any of its main conclusions. 
It is likely that some of these choices have even better properties than the ones that 
were taken. Perhaps a relatively simple change could be used to produce an effective 
metric gab that satisfies the vacuum Einstein equation beyond the regime of linearized 
perturbation theory. 



Appendix A. Generalized Killing fields 

The notion of a generalized Killing field (GKF) used in this paper was developed in 
pS] . where such objects were referred to as Killing- type generalized affine coUineations. 
Their main properties are summarized in Sect. |2.2[ For completeness, this appendix 
provides explicit definitions. It is heavily based on the description in |40j . 

Everything in this appendix is formulated on a spacetime {J^,gab)- Note, 
however, that the main text discusses GKFs constructed using different metrics. 
Besides the geometry, a generalized Killing field Sq also requires for its construction 
a smooth timelike worldline F = {7s |s G M} and a future-directed timelike vector field 
n° G Tj^-M defined along F. Note that s is not required to be proper time and 
needn't lie tangent to F. 

A specific generalized Killing field may be fixed by choosing a time sq together 
with tensors Aaiso) and Bab = S[afc]('So) at 7^0. The Killing transport equations 

^Aa{s)-j%a{s)^0 (A.l) 

as 

^Babis) + Rab/{ls)rsAd{s) = (A.2) 
as 

are used to uniquely extend these tensors to all of F. Note that the skew symmetry 
of Bab is preserved by this prescriptiorfpf] 



ft It is possible to use initial data for which Bab is not skew. The vector field that eventually results 
generalizes a homothety or other non-Killing afiine coUineation j25j . 
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Now consider all pairs (7s, w"), where v"" € T^^M is orthogonal to n^. This forms 
a subset T±T of the tangent bundle TAi. For any element of T±r, one may associate 
an afhnely-parameterized geodesic y{w) whose initial point is y{0) = 7^ and whose 
initial tangent is y"(0) = As long as these geodesies can be extended sufficiently 
far, the map (7s, — y{l) is a smooth function from T±T to A4. Its Jacobian is 
clearly invertible at (least at) every point (7^, 0), so it follows from the inverse function 
theorem that the given map defines a diffeoniorphism on some neighborhood W of F. 
This will be the region in which the GKFs are to be defined. It is assumed in the 
main text that the body whose motion is being studied always lies inside this region: 
W C W. 

We now define the GKF Ea{x) associated with a choice of Aa{so) and Bab{so)- 
The diffeomorphism just described may be used to uniquely associate x with some 
{IstI'"') € T±T. Use this pair to construct a geodesic y{w) as before. The GKF is then 
be computed along y{w) by solving the Jacobi (or geodesic deviation) equation 

^ - Rab/y'r^d = 0, (A.3) 

with initial data 

Sa(7«) -A(s), (A.4) 

^^^=v'Bba{s). (A.5) 
aw 

The given equations uniquely define Sa(a;) throughout W once Aa and Bab are given 
at any one point on F. More detailed discussions may be found in [TUlfSSlUn] . 
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